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Miles above the earth, the slender “Wac Corporal” rocket 
blasts loose from a modified German V-2, and attains a 
speed of 5,000 mph and an altitude of 250 miles—fastest 
and highest flight by any rocket. This most successful 
firing, from White Sands Proving Ground, New Mexico, 
February 25, 1949, gave scientists important information 
on little known phenomena of outer space and other data 
required in missile research. 

The “bumper” rocket used in this spectacular test was 
developed by Army Ordnance working with General 
Electric, California Institute of Technology Ballistics 
Research Laboratories, and Douglas Aircraft Company. 
Structural design of the “bumper” and construction of the 
“Wac” were among the contributions made by Douglas. 

This is just one of the many missiles projects, so vital 
to America’s defense, to which Douglas engineers have been 
devoting their skill and energies since 1941. 


DOUGLAS AIRCRAFT COMPANY, INC., SANTA MONICA, CALIFORNIA 
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GETTING YOUR ANSWERS 


IBM’s vast engineering know-how is helping American business, 
industry and the Armed Forces get the answers . . . fast. Through 
its leadership in applying electronic principles to calculators and 
other types of punched card business machines, IBM has given 
greater speed, accuracy and economy to the nation’s vital processes 
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of calculating and accounting. 


Already thousands of IBM Electronic Business Machines are in 


everyday use. We are continuing to manufacture them in quantity 


... as fast as quality production will permit. 
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Optimum Thrust Programming for a Sounding Rocket 


By H. S. TSIEN! and ROBERT C. EVANS? 


Tie problem of optimum thrust programming for a 
sour ling rocket of minimum starting weight to reach 
spec fied height with given final weight and propellant 
char :cteristics is first formulated as a problem in varia- 
tions! calculus. The general solution for arbitrary drag 
func ion is given. The solution is then applied to two 
special cases, one with quadratic drag dependence on 
velocity and the other with linear drag dependence 
on velocity. Complete numerical data are given. The re- 
sults are then compared with the results of constant 
thrust to show the advantages of thrust programming. 
Thrust programming is shown to be able to increase ap- 
preciably the pay load of a high altitude sounding rocket. 


KF )R a rocket in vertical flight, the aerodynamic 
drag of the rocket body and the gravitational pull 
are in the same direction and directly opposite to the 
direction of the thrust force. The performance compu- 
tation is thus relatively simple. It was indeed carried 
out by F. J. Malina and A. M. O. Smith (1). They 
have shown that if the aerodynamic drag is absent, then 
the best way of using the propellant is to use it in the 
shortest possible time. Theoretically, a rocket of given 
weight fraction of propellant will reach the highest alti- 
tude if the thrust is applied as a single impulse and thus 
reaches maximum velocity immediately. This result 
can be understood, perhaps, by considering the other 
extreme of having the thrust equal to the weight of the 
rocket at every time instant. Then the rocket, having 
zero acceleration at all times, will not be able to leave 
the ground. Long drawn-out operation of the rocket, 
therefore, is definitely unprofitable if aerodynamic drag 
can be neglected. 

When there is aerodynamic drag, the high velocity 
reached at low altitudes by an impulsive start will give 
avery high drag which tends to reduce the maximum 
altitude of the rocket. In fact, the calculations made 
by Malina and Smith, and more recently by Ivey, 
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Bowen, and Oborny (2), show that, assuming constant 
thrust, the optimum initial acceleration is one to three 
g, depending on the ratio of drag and weight of the 
rocket. The more general problem, however, is that of 
optimum thrust programming, i.e., to find the optimum 
time variation of the thrust for maximum altitude. This 
is the theoretical optimum design. Practical design is 
complicated by the added weight on the power plant 
to make the thrust variable and certainly cannot reach 
the theoretical optimum condition. The theoretical 
optimum gives a standard of comparison and shows 
how much could be expected by varying the thrust. 

The problem of optimum thrust programming was 
studied quite early by G. Hamel (3). He made the 
simplifying assumptions that the density of air de- 
creases exponentially with altitude, and that the ef- 
fective exhaust velocity c of the rocket motor and the 
gravitational constant g do not vary with altitude. 
Using variational calculus he gave the solution of the 
problem of optimum thrust programming. However, 
his paper is very brief and is not easy to understand. 
It is the purpose of the present authors to give a com- 
plete discussion of the problem together with more ex- 
tensive numerical data. 


Formulation of the Problem 


Let M be the mass of the rocket and s the height at 
time instant ¢. Following Hamel, the effective exhaust 
velocity c of the rocket motor will be taken as a constant. 
The velocity of the rocket is ds/dt and will be denoted as 
8. Foraspecified rocket body, the aerodynamic drag D 
is a function of the altitude s, and the velocity s. If g 
is the constant value of the gravitational force per unit 
mass, the equation of motion of the rocket during the 
powered flight is 


dM , M \ds 
dt a c 


The altitude s is measured from the point of launching. 
Therefore, at the beginning of powered flight, ¢ = 0, 
s = 0, M = M°,and 8 = 8) where M° is the mass and & is 
the velocity after boosting. At the end of the powered 
flight, = M = M,, s = = 3. is then the 
final mass. Equation [1] together with these boundary 
conditions gives the following expression for mass M°, 
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M® = exp (- *) x 


0 


If M, is the initial mass of the rocket including the 
propellant for boosting by a single impulse to the veloc- 
ity so, then 


Therefore, by combining [2] and [3], 


Equation [4] does not explicitly contain the velocity 
after boosting, or initial velocity %. If the initial 
velocity is zero, the acceleration of the rocket is smooth. 
If the initial velocity is not zero, the rocket starts with 
an impulse and then gradually accelerates. In any case, 
however, Equation [4] gives the total initial mass in- 
cluding the boosting charge. 

Now let the problem of optimum thrust programming 
be formulated as follows: Given M,, c, g, and the drag 
function D(s, 3), what should be the function s(¢), such 
that M, is a minimum? The auxiliary conditions are 
s(0) = O and that s, and 8, must be such as to reach the 
specified summit altitude S. To reach the given sum- 
mit altitude with the specified M, and D(s, 8), s; and 3 
are related, say 


where ¢ is a given function. For instance, at very high 
altitudes where the aerodynamic drag is negligible due 
to the low air density, 


[6] 


To find the conditions for the solution of this varia- 
tional problem, let the required s function be 


with s(0) = 0. 
Now let there be an arbitrary function y(t) such that 


but otherwise completely unspecified. Then the “neigh- 
boring” functions to s(t) can be constructed as 


S(t) = a(t) + 


where k is a parameter but not a function of time. Be- 
cause of Equation [8], § satisfies the initial condition 
5(0) = 0. The duration of powered flight, or the burn- 
ing time, for the optimum solution is t;. For the neigh- 
boring solution, the burning time is 4; + «. Thus k isa 
function of «. For the optimum solution, k and ¢ both 
vanish. Therefore, 


k(0) =0 


By considering only terms up to first order in e, 
= Ht + €) = + € &(t) + n(th)..... {11] 
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However and §, must satisfy Equa‘ ion 
[5] so that the neighboring solutions will repre: ent 
rockets reaching the specified summit altitude, the 


stated auxiliary condition. Therefore, 


where = 


3, and §, from Equations [11] and [12] can be sulsti- 
tuted into Equation [13]. After some simplifica' ion, 


one has 


This equation then determines the value of k’(0). \ith 
the k’(0) so determined, the neighboring solutions will 
definitely satisfy all the auxiliary conditions. 

With »(t) specified, the total initial mass Mo wii! be 
dependent upon e. Let 


= Ble, (=) [15] 
Then by substituting Equations [9] and [12] into |4], 


ite 
= F + + Rea, 0) dt + 


(th) + € (th) + + gh + (16] 


M, exp 


The condition for s(t) to correspond to the optimum 
solution can now be stated simply as 


7 


By carrying out the required differentiation, one has 


oF 
k (0) n(ti) + F(s:, $1, th) + 


+ gl 
But »(¢) is arbitrary other than the condition (0) = 0. 
Therefore, in order for the above expression to vanish, 


oF d {oF 
ds dt [18} 


and 


k/(O)n(ts) (55), + + 


My [& + 9 + exp =0....[19) 


Equation [18] is the familiar Euler-Lagrange differential 
equation. Equation [19] is the result of the auxiliary 
condition of the problem. 

By eliminating k’(0) between Equations [14] and 
[19], one has 
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and 


4, 4) + Mi [} (3). — (4 +g) + 


But 7(/) is arbitrary. Therefore, in order for the above 
equation to be true, the sums of quantities which 
multiply into n(t) and 4(t,) must be zero separately. 
There! ore 


M, (7). (8, + g) exp (ot = 0.... [20] 


and 


F(s, 4) + (3 exp (+. + = @....[21] 


Equa ions [18], [20], and [21] are now the complete 
answir to the variational problem. That is, the thrust 
progr imming must be such that Equation [18] is satis- 
fied «i every instant of the powered flight and, in ad- 
dition, at the end of the powered flight the conditions 
given by Equations [20] and [21] must be fulfilled. 
These conditions can be reduced to simpler forms if 
the relation in Equation [15] is reintroduced. Then 


Equation [18] becomes 

oD oD. 

1joD 


2 (8 + . - [22] 


(2k + 9) +2 


When the drag is specified as a function of s and 
Equation [22] gives the differential equation for the 
trajectory s(t). Equations [20] and [21] now become 


(52), + +9) = 
and 
D(s:, + [ (72). & + [24] 


where the subscript ( ); denotes the quantities evaluated 
att = t. Equation [24] is, however, automatically 
satisfied if Equation [5] represents the relation at the 
beginning of coasting flight. The reason is as follows: 
During coasting, the rocket motor is stopped, and there 
isno expenditure of propellant; therefore, dM/dt = 
Then Equation [1] reduces to 


This is the same as Equation [24]. 
By eliminating (“*) between Equations [23] and 


24], the condition at the end of the powered flight is 
finally expressed as 


+ Di + Mg) \(2 
+ + Mig) (Di + Mig) 
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where D, = D(si, 8:) is the drag at the end of burning. 
However, the factor (M,8; + D,; + Mg) is never zero, 
therefore, 


+ = D; + [26] 


The problem of optimum thrust programming can now 
be discussed in more concrete terms. Since the aerody- 
namic drag D enters linearly and homogeneously into 
Equation [22], that equation is actually a second-order 
differential equation for s(t), independent of the size of 
the rocket body. However, being a second-order dif- 
ferential equation with only one initial condition s(0) 
= 0, the initial velocity or velocity after boosting is 
yet free and undetermined. It is determined, however, 
by the condition, Equation [26], at the end of the 
powered flight. In other words, for a given size of the 
rocket and for a given final mass M,, there is a corre- 
sponding optimum booster velocity 8 and subsequent 
optimum thrust programming for any specified summit 
altitude S. In general then, the optimum solution al- 
ways involves an impulsive start from rest. -This is a 
characteristic of the problem. 


Quadratic Drag Law 


To carry out the computation, the air density is as- 
sumed to decrease exponentially with respect to alti- 
tude. As a first example, the drag is taken to vary as 
the square of the velocity. Then the aerodynamic drag 
is given by 

= exp [27] 


This corresponds to a constant drag coefficient. Sub- 


stituting this expression into Equation [22] yields 


where 
B [29] 
both are nondimensional parameters. Integrating 


Equation [28] for ¢(v) and s(v), 


+(1 — Bw — 26 
2 + (1 — — 
as =v—v + 
where : 


and In indicates the natural logarithm to the base e. 
By using Equations [27], [28], [30], and [31], the 
mass M at any time instant is found to be as follows: 
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2 
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+ (1 — — 26 = 
+2 


Pig v = uv, t = 0 gives the mass M° after boosting: 


vo? + (1 — — 28 
Uo + + 2) vn? + (1 — Bn — “A 


exp (» ... [34] 


Mo 


The initial mass of the rocket including the booster 
change is given by Equation [3], or in terms of the non- 
dimensional velocity v, 

My = M° exp (vo) 


Therefore, 
[35} 
The thrust at any instant can be obtained from 
Equation [1] noting that F = c 
ig exp (—as) + (1+ *) [36] 
where the acceleration §/g can be computed from Equa- 


tion [28]. 

The end conditions are given by Equations [5] and 
[26]. Substituting Equation [27] into Equation [26] 
gives one condition 


This equation can be easily integrated. By using the 
condition that at the summit, s = S, § = 0, one has the 
following equation for the velocity v; at the beginning of 
the coasting flight 


2 = — 28 exp (28 exp (- 28 ve .139) 


where 


& = exp (—as;), é& = exp(—aS)........ [40] 


The integral in Equation [39] can be evaluated using 
the series expansion 


ff exp (- az) @ 


With any fixed value of S and the drag parameter 
Woc?/Mig, Equations [37] and [39] determine and 
as;. The parameter Woc?/M,g is the nondimensional 
drag and weight ratio of the rocket. 

Calculations were carried out for two sets of summit 
altitudes and exhaust velocities assuming a = 1/22,000 
ft. One case was for a summit altitude of 500,000 ft and 
an exhaust velocity of 5500 fps. The other case was for 
a summit altitude of 3,000,000 ft and an exhaust ve- 


(ax)? 


+5 
.. [41] 


=In|z|-j 


locity of 8000 fps. Actually, an iteration procedure was 
adopted to fit the end conditions. For any chosen 
value of the drag weight ratio Woc?/Myg, a value of », 
was first assumed and substituted in Equation [37! to 
solve for as;. This result was put into Equation [°9] 
which was solved for 1;._ This process was repeated u: til 
the desired accuracy was obtained. A plot of Woc?/.\/g 
versus v;, which was found to be almost a straight | ne 
on semilogarithmic graph paper, with v; plotted on he 
linear scale, simplified the process by giving an excell:-nt 
first estimate of the velocity. The initial velocity rz tio 
vo was then solved for bya trial and error procedure us ng 
Equation [31]. 

Values of Woc?/M,g were assumed for the calculatic as, 
but the final results were put in terms of My. This ‘as 
done because physically the results in terms of he 
initial mass are easier to visualize. The quan: ity 

‘yc? is the drag of the rocket body at a flight veloc ity 
equal to the exhaust velocity c and at sea level. “ ‘he 
parameter then expresses the ratio of this drag to ini ial 
weight. The results are shown in Figs. 1 to8. The «iis- 
cussion of these results is postponed until a later :ec- 
tion. 


Linear Drag Law 


For supersonic flight of a rocket, perhaps, a better 
approximation of the drag is given by 


= (Aé + B) exp (—as)............. [42] 


This again corresponds to an exponential air-density 

law, but the drag coefficient decreases as the velocity 

increases. Substituting Equation [42] into Equation 

[22] gives 

_ + (B/Ac — By — [8 + B/AC1 + 
A(v + 2 + B/Ac) 


where v and £ are the same as the previous section and 
are defined by Equation [29]. Equation [43] can be 
integrated and gives 


ws v? + o(B/Ac — B) — [68 + B/Ac(1 + 8B)] 
vo? + vo(B/Ac — B) — [8B + B/Ac(1 + B)] 
68+ B/Ac) 2v + (B/Ac — B) — 
2r 2v + (B/Ac — B) +A 
_ 2v0 + (B/Ac — B) +d [4] 

+ (B/Ac — 8) 


and 
as =v—%m+ 


(1 4) u(B/Ac — 8) — [B + B/Ac(1 + 8)] 
2 vo? + vo(B/Ac — B) — [8 + B/Ac(1 + 8)] 

B+ B/Ac) + (B/Ac — 6) — 

2v + (B/Ac — B) 
+ (B/Ac B) 145) 


where 


h = V(B/Ac + B)* + 4(B/Ac + B)........ | 46] 


The mass M at any instant is found to be 
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44] 


45] 


46] 


a -B(exp vo) + vo(B/Ac — B) — (B + B/Ac{1 + 8])} X 


v+ B/Ac+1 
jo? + o(B/Ac — B) — (8 + B/Ac{1 + 
vy + B/Ac+1 


n® + 0(B/Ac — 8) — (8 + [1 + * 
exp ¢ + exp + 


The initial mass My including the boosting can be ob- 
tain d from the above equation and Equation [35] 
= exp + + 

8 [tm + B/Ac +1) + B/Ac +1) X 


+ v( B/Ac — 8) — (8 + B/Ac{1 + 
— »(B/Ac — B) — (8B + B/Ac{1 + 8B}) 


-exp (vo). . . [48] 


The thrust can be obtained from Equation [1]. Thus 


ir (1 + [49] 


hg 
The xcceleration s/g can be determined from Equation 
[43 | 
The end condition is given by Equations [26] and 
[5]. Substituting Equations [42] and [43] into Equa- 
tion |26], one has 


fw + B/Ac) exp (—as) + 


a + B/Ac — B/Ac) exp (—as) = 1... [50] 


To supply the relation indicated by Equation [5], the 
coasting flight trajectory must be determined. By 
using Equation [42], the equation of motion for the 
coasting flight is 


Ac 
¢ + 1) + Ma (v + B/Ac) exp (—as) = 0..... [51] 


This equation cannot be solved by simple quadrature. 
However, the effect of air drag is generally small. Then 
Equation [51] can be solved approximately by first 
neglecting the air drag and then making the necessary 
small correction. Let v® be the value of the velocity 
ratio without air drag and v’ be the small correction on 
Then 


The equation for v® is simply 


d(vo)? 29 
ds c 


Now substitute Equation [52] into Equation [51] and 
retain oo linear terms in v’, 


0 + B/Ac) exp (—as) = 0... [54] 


This is the differential equation for v’, the correction 
term. Equation [53] has the solution 

= V2B(aS — [55] 
Substitute Equation [55] into Equation [54], 


dg 2¢ Mig ad exp (¢ aS) 
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where 


¢ = aS — as 
Solving for v’ 


To obtain a suitable series for the computation of 
the integral in Equation [56], consider 


05) = f° exp (oar 
Let x = ¢(1 — w), then 
vi — wexp (— fu) du 


“texp (x) dz . [56] 


f(g) exp 


By expanding the radical in the integrand into a power 
series and integrating term by term, f(¢) can be shown to 
be equal to the following series 


f(t) = Vi [ - fexp (¢) - 


1 
Vig MG) + B/Ac 


(exp DEX 
exp (— a@S).. [57] 


1 
V28 


Equations [50] and [57] must be solved simultaneously 
to determine the velocity and altitude at the end of 
burning for a given drag weight ratio, Ac/M,g, for the 
rocket. The quantity Ac is now the drag of the rocket 
body at a flight velocity equal to the exhaust velocity c 
and at sea level, if B = 

Calculations were carried out for a summit altitude of 
3,000,000 ft, an exhaust velocity of 8000 fps, and no 
drag at zero velocity, assuming a = 1/22,000 ft. An 
iteration procedure identical to that used in the pre- 
vious example was used. The results are shown in Figs. 
9 to 12. The parameter here is Ac/Mog, the ratio of 
the drag of the rocket at velocity c and at sea level to 
the initial weight. 

In order to determine the gain possible by using the 
optimum thrust programming as against the conven- 
tional constant thrust rocket, a few cases of constant 
thrust were computed using the same drag laws. 
It was assumed the constant thrust rocket had no 
initial velocity and the ratio of thrust to initial 
weight F/Mog was 2.7. This was chosen to agree 
with the best initial acceleration given by Ivey, 
Bowen, and Oborny (2) for high-performance rock- 
ets. The comparisons were based on equal ratios 
of the drag at the exhaust velocity c and at sea level to 
the initial weight Mog and summit altitude S. The mass 
ratio My/M, is now larger than the optimum value due 
to the fact that the final mass M, is now smaller. A 
step-by-step procedure was used in computing the mass 
ratio of final mass to initial mass which gave an error of 
approximately one per cent. This error was always 
nonconservative and would therefore favor the con- 
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stant-thrust examples in the comparisons. The results 
are indicated by the dots in Figs. 1, 2, and 9. 


Discussion of Results 


Both the linear drag law and the quadratic drag law 
cases have certain common characteristics. As seen 
from Figs. 3, 4, 7, 8, 10, and 12, for low values of the 
drag parameter, the average thrust and initial velocity 
were high. Therefore, the rocket will burn its fuel 
rapidly and less energy will be expended against gravity. 
As the drag parameter increases, more energy is used to 
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overcome the drag, and it becomes advantageous to 
have a lower velocity at low altitudes. Therefore, the 
initial velocity and thrust will be lower, and the altitude 
at the end of burning will be greater than in the casv of 
low drag parameter. But it is important to note that 
there is always a finite initial velocity so, so that the 
optimum thrust programming involves a boosting im- 
pulse. Another result the two drag law cases have in 
common is the effect of the drag parameters on the 
ratio of the initial mass to the final mass, M)/M, (Figs. 
1, 2, and 9). As the drag parameters increase, ‘his 
ratio M,/M, increases for a given summit altitude. As 
the drag parameters decrease, the ratio M)/M, asyinp- 
totically approaches the no-drag value of exp so/¢ where 
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is such that the rocket has the kinetic energy at zero 
altitude equal to the potential energy at the summit alti- 
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In general, for optimum programmiag, the thrust 
after the initial impulse should increase with altitude. 
Exceptions to this occur for portions of the powered 
flight when the values of the drag parameter are very 
high, or when the rocket has a very high performance. 
In these cases, for a portion of the burning time, the 


thrust can be seen to decrease. This is the result of 
either a very slowly increasing acceleration of the 
rocket or a very rapidly decreasing mass of the rocket. 
Since the accelerating force is the major fraction of the 
thrust and equal to the product of the instantaneous 
mass and acceleration, the thrust can decrease. In 
general, however, the thrust should increase. For a 


constant geometry rocket motor, this is accomplished 
to a limited extent by the decreasing atmospheric pres- 
sure acting on the nozzle exit. But for most cases, such 
an increase is only a fraction of the ideal, as determined 
by the present calculation. 
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The linear drag law case has a higher initial velocity 
and shorter burning time than has the quadratic drag 
law case for similar mass ratios My/M;. This is true 
because the drag does not increase as rapidly with 
velocity in the linear drag law case; therefore, the 
linear drag law case will favor the higher velocities and 
shorter burning times. The drag in the quadratic law 
case increases as the square of the velocity; therefore, 
this case will favor the lower velocities and longer burn- 
ing times. 

To give an example of the magnitudes of the ratios 
Woc?/Mog and Ac/Mvg for the drag and the initial 
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weight, consider the case of the V-2 rocket. At a veloc- a 
ity of 6000 fps the drag is roughly 100,000 lb. The tl 
initial weight of the V-2 is approximately 25,000 lb. a 
This gives a ratio of drag at the exhaust velocity at sea tl 
level to the initial weight of roughly 4 if the exhaust nr 
velocity is 6000 fps. Therefore, 4 might be an approxi- ir 
mate value of the ratio Ac/Mog. In the quadratic crag 

law case, a better estimate of Woc?/Mog might be 10, fo 
This will account for the quadratic drag law giving a su 
lower value of the drag than is actually the case during lo 
the major portion of the powered flight, especially at the F “4 
lower altitudes. Since the summit altitude of the V-2is J 
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approximately 100 miles, Fig. 3 can be used to estimate 
the optimum booster velocity to be roughly 1500 fps, 
assuming Woc?/Mog = 10. If the drag is relatively less, 
the booster velocity will increase. Therefore, the opti- 
mum thrust programming involves considerable boost- 
ing. 

Several constant-thrust mass ratios were calculated 
for comparison with the ideal mass ratios. It was as- 
sumed the constant thrust rockets had no initial ve- 
locity and followed the same drag laws as the ideal 
cases. An initial acceleration of 1.7g was chosen to 
agree with the best initial acceleration given by Ivey, 
Bowen, and Oborny (2) for high performance rockets. 
This value might not be the optimum for the constant 
thrust cases computed because it is dependent on the 
drag laws followed by the rocket. The values for the 
constant thrust rocket calculations can be seen in Figs. 
1,2,and9. The advantages of the ideal rocket would be 
less if the constant thrust rocket had the best possible 
initial velocity and acceleration. In fact, it seems cer- 
tain that the initial acceleration of 1.7g is too small for 
the linear drag case of Fig. 9. The possible gain by 
using optimum thrust programming is thus more 
correctly represented by Fig. 2. Then it is seen that 
the gain is small in terms of the initial mass Mp. On 
the other hand, one must bear in mind that the pay load 
of a high altitude sounding rocket is only a very small 
fraction of the initial weight. Therefore, even if the 


gain by better thrust programming may be small in 
terms of the initial weight, the gain in terms of pay load 
could be appreciable. Of course, the numerical results 
obtained by the present calculation are only qualitative. 
For any practical design, besides consideration for the 
additional weight needed for thrust variation, the 
optimum-thrust programming has to be determined 
with a more factual drag versus velocity relation. 
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A Theory of Unstable Combustion in Liguid Propellant 
Rocket Systems 


By MARTIN SUMMERFIELD 


Member ARS, General Editor, Aeronautics Publication Program, Princeton University, Princeton, N. J. 


On the basis of an hypothesis that low-frequency oscilla- 
tions (chugging) sometimes observed in liquid propellant 
rocket engines, are the result of oscillatory propellant flow 
induced by a combustion time lag, conditions for the sup- 
pression of such oscillations are derived. It is found that 
stability can be achieved by increases in the length of feed 
line, the velocity of the propellant in the feed line, the 
ratio of feed pressure to chamber pressure, and the ratio 
of chamber volume to nozzle area. Equations are given 
for the frequencies of oscillation. Examination of the 
equation for stability indicates that self-igniting propel- 
lant combinations are likely to be more stable than non- 
self-igniting systems. 


Nomenclature 


Physical Quantities 


A = cross-sectional area (normal to the flow direction) 
c* = characteristic velocity of the propellant 
K = over-all frictional pressure loss coefficient 


1 = length of flow passage 

L* = characteristic length of the rocket motor 

mass flow (mean value) through the rocket system 
p = pressure (absolute) 

Ap = pi — pe = pressure difference from tank to chamber 
R = gas constant of combustion gas (based on unit mass) 
T = temperature (absolute) 

if 

B 

w 

p 


= 
=. 
= 


= velocity of flow 

= volume of chamber 

= volume compressibility of liquid 

= damping factor of oscillation 

= frequency of flow oscillation (radians/sec ) 
= density of the liquid 

= combustion time lag 


Subscripts 


, refers to propellant tank 

2 refers to feed line 

3 refers to injector orifice 

- refers to combustion chamber 

: refers to exhaust nozzle throat 

refers to equilibrium value of v 

9 refers to combustion gas 

Parameters 

A, B,.C, a, b, c, U, u, x, 6 are employed in the analysis to simplify 
algebraic manipulation and are defined in the text. 


Introduction 


IQUID propellant rocket engines of all types are 
generally designed by the engineer to deliver a 
fairly constant steady thrust for a duration that may 
extend from several seconds to as much as an hour. 


Received April 16, 1951. 
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However, in many cases it has been reported that the 
desired steady operation does not occur in actual test, 
and instead a condition variously described as “‘rough 
burning,” “‘chugging,” “screaming,” or simply unstzble 
combustion may take place. Frequencies ranging 
from 10 cycles per see to as much as 5000 cycles per sec 
have been observed in oscillographic chamber pressure 
traces, amplitudes from a few per cent to as much as 
50 per cent of the mean chamber pressure, and in many 
cases, the oscillation was not truly periodic but seeined 
to be merely a series of random fluctuations. Not 
only is the resulting thrust vibration undesirable from 
the standpoint of possible damage to the structural 
elements or instruments in the vehicle, but in extreme 
cases, failure of the power plant itself can occur (1).! 

The theory presented here does not attempt to ex- 
plain all cases of unstable combustion, but it may pro- 
vide an explanation for one type of instability, perhaps 
the most serious type, sometimes called chugging. 

The phenomenon of chugging was first observed by 
the author in October, 1941, during a series of tests of 
a 1000-lb-thrust, nitric-acid—gasoline rocket motor at 
the Jet Propulsion Laboratory of the California Insti- 
tute of Technology. Upon ignition, combustion would 
proceed smoothly at first, but would rapidly become 
rough, and in many instances, after five or ten seconds, 
severe fluctuations in chamber pressure would be 
taking place. Frequently, the bolts of the combustion 
chamber would rupture before the propellant flow could 
be turned off. 

An explanation for this oscillation was advanced on 
the basis of an hypothesis that a time lag existed be- 
tween an arbitrary fluctuation in propellant flow and its 
subsequent manifestation in the combustion chamber 
pressure. This time lag is the result of the finite rate 
of the over-all combustion process, and is determined, 
therefore, by the kinetic rates of mixing, vaporization, 
and chemical reaction. The argument was pursued to 
show that the instability could be suppressed by increas- 
ing the pressure drop across the injector. A numerical 
illustration of this argument may be helpful for an 
appreciation of the analysis below. 

A particular rocket engine with a compressed-gas 
type of feed system is designed to operate with a cham- 
ber pressure of 300 psia.. In the first case, assume that 
the feed pressure is 500 psia. Now suppose a momel- 
tary decrease in chamber pressure occurs after steady 


! Numbers in parentheses refer to references on page 114. 
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operation is achieved, the decrease being from 300 to 
200 psia. Assuming a square-law pressure drop across 
the injector, and neglecting the inertia of the liquid in 
the feed line, the flow rate will increase by 23 per cent 
above the design value and will subsequently (after 
the time lag mentioned above) produce a 23 per cent 
increase in chamber pressure above the design value, 
namely, 369 psia. The injector pressure drop is now 
only 131 psia, the corresponding flow rate is therefore 
reduced to 81 per cent of the design value, and the 
chaniber pressure will then fall (after the time lag), to 
81 per cent or 243 psia. The flow rate then increases 
again, and it is possible to continue the calculation 
inthe same manner. In this case, after one cycle, it is 
apparent that the disturbance is decaying, the ampli- 
tude having decreased from 100 psia to 69 psia to 57 
psia. 

Next, consider a second case in which the design 
chamber pressure is also 300 psia, but the feed pres- 
sure is only 400 psia.- If an arbitrary decrease of 100 
psia is assumed as above, the sequence of chamber 
pressure values becomes 300, 200, 420, 0. In this case 
it is clear that the fluctuations will increase in strength 
as time proceeds. A third case, assuming a feed 
pressure of 450 psia, exhibits neutral stability, i.e., 
the amplitude of the fluctuation remains approximately 
constant. 

It appears from such crude considerations that the 
injector pressure drop is a controlling parameter. — It 
will be seen below that this numerical conclusion can 
be generalized in the following statement: Instability 
isnot possible if the injector pressure drop exceeds one 
half the mean chamber pressure. However, it will 
be seen that the converse is not always true. That is, 
instability can be suppressed even when the injector 
drop is less than half the chamber pressure. 

Of course, even in its qualitative form, the theory 
suggests that instability can be eliminated by reducing 
the time lag in the combustion process. This was 
accomplished, in the 1941 program, by switching to 
aniline as a fuel instead of gasoline, and indeed, with 
no change in design or operating conditions, smooth 
running was achieved. 

In formulating the problem in 1941, the author’s 
analysis neglected the inertia of the liquid propellant 
in the feed lines.2 The flexible coupling responsible 
for this self-excited oscillation was assumed to be the 
tapacitance of the combustion chamber. A_ recent 
paper by Gunder and Friant (2) presented an analysis 
of chugging in which the essential hypothesis was also 
the time lag described above (it was independently 
tonceived by them), but instead of the chamber 
tapacitance the important term was the inertia of the 

*The author is indebted to Dr. Theodore von Karman for 
profitable discussion that led to the basic hypothesis of a com- 
bustion time lag and for suggestions on the mathematical ap- 
ot to the problem. (Private communication, November, 
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liquid in the feed lines. It is shown in their paper that 
instability is not possible when the pressure drop exceeds 
half the chamber pressure. Yachter and Waldinger, 
in an unpublished communication, 1949, considered 
the case in which liquid inertia is neglected, and reached 
the same conclusion. 

Since similar conclusions resulted either from a con- 
sideration of liquid inertia alone or chamber capacitance 
alone, it appeared logical to carry out the analysis with 
both effects present. This is done below, and not only 
does it include each of the above analyses as special 
cases, but certain conditions for stability can be derived 
in such simple form that physical conclusions suggest 
themselves immediately. 

In passing, a few remarks concerning other types of 
nonsteady operation are pertinent in order to clarify 
the assumptions of the present theory. It is assumed 
here that the pressure is uniform throughout the com- 
bustion chamber at any instant, that is, that the period 
of oscillation is long compared with the time required 
for a pressure wave to traverse the combustion chamber 
in any direction. However, high-frequency screaming 
has been observed in some rocket motor tests, which 
indicates that chamber gas oscillations are possible, in 
addition to the system oscillations treated here. These 
high frequencies usually agree fairly well with those 
expected on the basis of axial resonance (organ pipe 
notes) or transverse resonance. What is not clear, 
however, is the nature of the process that energizes 
such oscillations and why they appear only under 
certain conditions of injection or chamber configura- 
tion. It is suspected that the oscillation is energized 
by a proper interaction between the variations in pres- 
sure and the rate of combustion (that is, the time lag 
is some function of the instantaneous pressure), but 
the exact formulation of this interaction is not known.’ 
The remedy is not clear unless it is to try arbitrary 
changes in the spatial combustion pattern. These 
high-frequency vibrations are usually of small ampli- 
tude but are generally accompanied by significant 
increases in heat transfer to the chamber walls. The 
danger of this type of oscillation lies in the possibility 
of a burn-out of the chamber. 

A third category of unsteady operation is connected 
with either mechanical or hydraulic oscillations in the 
rocket system. Aside from the well-known _possi- 
bilities of structural vibration, valve flutter, feed pump 
pulsation, ete., there seems to be the additional pos- 
sibility of nonsteady hydraulic behavior of certain 
elements of the flow system, particularly the injector 
orifices. The danger of operating near the critical 
flow rate corresponding to the appearance of a vena- 
contracta in the injector orifice was first pointed out 
to the author by W. B. Powell‘ in 1942. The remedy 


3 An analysis based on this hypothesis has been advanced by 
Dr. L. Crocco of Princeton University. 

‘Jet. Propulsion Laboratory, California Institute of Tech- 
nology. 
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for this type of instability lies clearly within the domain 
of the hydraulic engineer. The present theory is not 
concerned with such effects. 


Analysis 


A schematic liquid propellant rocket system is shown 
in Fig. 1. It is assumed that the system is either a 
monopropellant type or, if a bi- or multipropellant 
system is under consideration, that the feed lines and 
injector orifices for the separate liquid reactants have 
identical hydraulic characteristics so that it behaves 
like a monopropellant system. Although it is possible 
to repeat the monopropellant analysis (below) for the 
case of a multipropellant system, the additional com- 
plication is hardly worth while until the basic ideas 
are confirmed by experimental checks. The physical 
principles and the resulting rules for overcoming 
instability will be qualitatively the same for the bi- 
propellant case. 

At any instant during the oscillation the rate of 
change of pressure in the combustion chamber is 
governed by the difference between the rate of evolu- 
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tion of combustion gas and the rate of gas flow through 
the nozzle.’ In accordance with the fundamental 
hypothesis of this theory, the rate of evolution of gas 
at a given instant is equal to the rate of flow of pro- 
pellant at an instant 7 sec earlier. 


dp. _ RT. RT. 
= — r) — L*c* [1] 


The effect of a changing chamber pressure on the 
flow of liquid propellant can be calculated by first 
setting up the general equation for the instantaneous 
flow rate with a specified difference between tank 
pressure and chamber pressure, and then differentiating 
this equation with respect to time. A convenient 
approach is to consider two control surfaces normal to 
the flow direction, one of area A; in the propellant 
tank upstream of the inlet to the feed line, and the 
other of area A; at the exit of the injector orifice, and 
then equate the total energy (work plus kinetic energy) 
entering the first control surface in a time interval 
dt to the energy leaving the second control surface 
plus the change of kinetic energy of the liquid contained 
between the two control surfaces plus the energy dissi- 
pated in the same time interval. 


(main + ; dt = +3 ps?) 
d (54 +5 Aolopv2” + 5 + 
(hit + + 
Carrying out the differentiation, replacing » and 


v3; by v through the continuity equation, the equation 
takes the following form: 


A? : - ]- 
— +500 [ (45 - - 


+h+h- . [2] 
As Ay 
In practical designs, < lo, and bs, 


Hence, without appreciable error, these terms may be 
neglected. The energy dissipations due to pipe line 
friction, valve losses, and imperfect flow in the injector 


orifices are contained in the term K 5 pv)”. 


Differentiating Equation [2] with respect to time, 
the desired equation is 
“pe + As? + k) + = 
Yquations [1] and [3] can be combined to eliminate 
p, and provide a single equation with a single dependen' 
variable, 7. Thus, differentiate Equation [1] to 
produce an equation containing the first and second 
derivative of p,. The first derivative can be replaced 
by Equation [3], and an expression for the second 
derivative is obtainable by differentiating Equation [3}. 
The equation for v2 is the following (dropping§ 


’ Symbols employed in the analysis are explained on page 108 
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subscript 2 of v2, and writing prime for the derivative 
with respect to time): 


RT ” A;? A;? , 
+ ape * K) vo!’ + k) (v’)? + 
+ K) + = 0. Al 


This nonlinear differential equation can be linearized 
by the following argument. If the amplitude of 
oscillation Av is small compared with the equilibrium 
value of v (which is a plausible assumption in testing 
for instability), then v, can replace v in Equation [4]. 
From Equation [2] 


p 


The same assumption that Av is small compared with 
y. makes it legitimate to throw out the (v’)? term in 
comparison with either the v’’ terms or the v’ terms. 
The latter comparison leads to this conclusion pro- 
vided RT./c*L* is identified as the maximum value 
of the frequency when conditions are near neutral 
stability. This identification appears later and pro- 
vides a posteriori justification for this approximation. 

Making the substitution, « = v’, the differential 
equation becomes 


In this equation 


To solve Equation [6], it is assumed that the general 
solution can be represented in the form of a linear sum 
of particular solutions, each of these being a sinusoidal 
oscillation modified by a damping factor.* 


n 


The successive values (A,, w,) are the roots of the 
characteristic equations obtained by inserting the 
above solution into Equation [6], and grouping real 
and imaginary terms: 

+ AX + B — w? + vos (wr) = 0...... [9] 
+ Aw — sin (wr) = 0..........[10] 


It is convenient to transform these equations by 
introducing new variables: 


wr = 
Ar 
[11] 
Br? = b 
Cr? =c 
Equations [9] and [10] become 
+ ar + b — 6? + ce-*cos6 =0......... {12] 


‘The method of solution employed here follows that of N. 
Minorsky for time lag problems (4). 
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(22 + — ce-7 sind =0............ 


By plotting these equations in the x, 6 plane, it is 
possible to locate all the intersections of Equations 
[12] and [13], and thus, for any given values of A, B, 
C’, r, to determine the behavior of the system. Those 
solutions involving positive values of x lead to in- 
stability; those solutions with negative x represent 
disturbances that decay. It is the objective of the 
designer to select values of A, B, and C such that all 
the terms in Equation [8] are damped (or at least the 
low-frequency terms, since the high frequencies are 
generally suppressed by viscous dissipation). 

A plot of Equations [12] and [13] would look like 
Fig. 2 for arbitrary values of a, 6, and c. Only the 
positive @ half of the plane need be examined. Both 
equations are invariant with respect to a change in the 
sign of 6. Equation [12] yields two branches, a positive 
one that is continuous and a negative one that exists 
only in certain ranges of @; Equation [13] produces a 
single-valued curve, also intermittent. The encircled 
intersections in Fig. 2 correspond to the terms in series 
[8]. The lowest frequency term is obviously the criti- 
cal one with regard to the over-all stability of the sys- 
tem. An unstable run may exhibit a finite number of 
low frequency modes, but the modes of higher frequency 
will be damped. 

Certain conditions can be derived for which solutions 
with positive values of x are impossible. Thus, if 
a > c, the curve given by Equation [13] would lie 
entirely in the negative x half of the plane. There- 
fore, this nequality is a sufficient condition for stability. 
The inequality can be expressed in terms of the ap- 
propriate design parameters by means of the defining 
Equations [7] and [11] and the following relationships: 


m = 


1 
2 ( 3? ) 


c* = p.A:/m 


Ap 


The condition for stability then becomes 


c*L* 2p 
RT. 


ly m 


Another useful form of this equation can be obtained 
by using the theoretical expression for c*: 


+20 

x 

| %.025 

sec. b=10.4 
a= 9.1 c=13.4 
12) 
-20 
FIG. 2. APPROXIMATE REPRESENTATION OF EQUATIONS [12] 
AND [13] 
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f(y) 1.55 for 1.1 << y < 1.3 


By substitution in Equation [15] the inequality takes 
the form 
m 
pe As 
Although this el expresses merely a suf- 
ficient condition for stability, it will appear in the follow- 
ing section, when the magnitudes of the terms are 
examined for typical cases, that the condition is suf- 
ficiently close to be useful as a design criterion. 
Another general condition for stability can be 
derived by combining Equations [12] and [13] and 
examining the possible roots (6,, 2») 


+ 6f(x) + g(x) = 0 \ 
f(z) = 2x? + 2ax + (a? — 2b) ......... [18] 
g(x) = (x? ax + b)? — 

The signs of f(x) and g(x) are of interest. The term 


(a2 — 2b) can be evaluated in terms of the defining 
Equations [7] and [11], and it can be seen that it is 
inherently positive. Consequently, in the positive 
x domain where our attention is focused, f(x) 
positive. Then, @? can be positive only if g(x) is 
negative, and conversely, it is possible to eliminate 
real, positive values of 6°, and hence of 6, by setting 
g(x) > O in the positive x domain. This condition 
can be assured by setting b > c. This inequality can 
be expressed in terms of the design parameters through 
Equations [7] and [11], yielding the following condi- 
tion for stability 


Pe 


As in the case of Equation [15], this inequality repre- 
sents merely a sufficient condition for stability. Clearly, 
stability is possible even if 2Ap/p, is less than unity, 
by reference to Equation [15]. 

The frequencies of the oscillations can be evaluated 
in principle by solving Equation [18] together with 
either of Equations [12] or [13]. However, considera- 
tion of Equation {18] alone can indicate the magnitude 
of the frequency of the particular mode in series 
[8] that makes its appearance as the stability 
conditions [15] and [19] are gradually relaxed. Thus, 
Equation [18] being quadratic in 6?, it can be solved 
explicitly for 6.,, the critically damped frequency 


ber = WV — 1/2f(0) + 1/2 — 49(0)........ [20] 


If it is supposed that the condition of critical damping 
occurs when (2Ap/p,) is near to unity, then by setting 
(c — b) < 1 and expanding the square root accordingly 
it can be shown that 

(i - 1/2 
[21] 


~ + 


As expected, w,, 0 as (2Ap/p.)er 1. 
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As a second case, it may be supposed that the condi- 
tion of critical damping is reached by relaxing condition 


[15]. Then, c may be replaced by a in Equation [20], 
and if it is further assumed that lm/p,As is of the 
2Ap 


aw ~ [ (282) (22° + 


By comparison with condition [15], it is apparent that 
in this case w,, is of the order of 277. 

A third case of interest may be developed by assum- 
ing that r < w,~!. This assumption makes it pos- 
sible to reconsider the basic differential Equation 6]. 
Thus, «(¢ — 7) can be expanded in a Taylor series: 


ut — r) = u(t) — ru’(t) +... 


, it can be shown that: 


same order as 


and the differential equation takes a well-known form 
+(A — Cry’ +(B4+Cyu = 0........ | 24] 


The condition for stability is that the coefficient 
(A — Cr) shall be positive; this provides the same 
relation as condition [15], but in this case it is a neces- 
sary condition, whereas condition [15] was merely a 
sufficient condition. The frequency in the critically 
damped situation is simply VB + C which reduces 
exactly to Equation [22]. (This result requires care- 
ful examination in each specific case, however, since 
VB + C is often of the order of 7~!.) 

By inserting typical design values in Equations 
{20}, [21], and [22], it may be concluded that frequen- 
cies greater than 200 cycles per sec are unlikely. ‘This 
supports the opinion expressed in the introduction, 
that this type of instability (chugging) is a different 
phenomenon from the high-frequency oscillations 
(screaming) sometimes observed. 

It is of interest to observe that the monopropellant 
analysis of reference(2) appears as a special case of this 
analysis. In Equation [6] allow L* + 0 


This equation can be identified exactly with Equation 
[5] of (2). Being of lower order, the differential 
equation for this limiting case can be solved for the 
critical condition in simple closed form (in contrast to 
the system of Equations [12], [13], [20] which do not 
yield so simple a solution); thus, it is unnecessary to 
resort to the relatively weaker sufficient conditions. 
The critical condition for stability then becomes 


lan (26 
4/1 — (2Ap/p.)? 


This is equivalent to Equation [13] of (2). This 
stability condition is useful whenever (lgr/p-Az) 
much greater than (c*L*/RT.)(2Ap/p.). It has 
the advantage of permitting the use of smaller values 
of (l.2/p-As) than the weaker condition [15]. 
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Discussion of Results 


Equations [15] and [19] lead to the following rules 
for overcoming instability in a liquid propellant rocket 
system: 


(a) Increase the pressure difference between supply 
tank and combustion chamber, either by reducing the 
area of the injector orifices or by inserting resistance 
elements in the feed lines. 

(}) Increase the volume or L* of the combustion 
chamber. 

(c) Increase the length of tubing from supply tank 
to combustion chamber. 

(/) Reduce the cross-sectional area of the propellant 
flow passages leading to the injector (or increase the 
mass flow per unit area). 

(-) In a bipropellant system, instability may occur 
in only one of the two feed systems. Therefore, an 
important step in eliminating instability is to increase 
the value of the left side of condition [15] pertaining 
to the weaker of the two. 

(/) Reduce the combustion time lag of the propel- 
lant. This may require changing over to a more 
reactive propellant or adding catalytic or combustion- 
promoting substances. This may be outside the 
bounds of permissible changes in practice. However, 
itis possible to change the configuration of the in- 
jector (somewhat empirically) or to preheat the pro- 
pellants to accomplish the same result. One way to 
accelerate combustion of the propellant is to provide 
a flow pattern in the combustion chamber in which 
hot gases are recirculated vigorously to mix with the 
incoming propellant. 

There is experimental evidence that some of these 
recommendations are in the right direction. It has 
been reported that an increase in Ap or an increase in 
L* has been effective in certain instances in eliminating 
instability (1). In the author’s experience, switching 
from nitric acid—gasoline to nitric acid—aniline in the 
identical rocket system removed the instability that 
had been encountered with the former. (See the 
remarks in the introduction.) In the case of nitro- 
methane, a redesign of the combustion chamber and 
relocation of the injectors so as to promote recircu- 
lation of hot gases and more prompt vaporization and 
reaction of the injected liquid has been found to sup- 
press instability. 

It is of interest to consider the numerical magnitudes 
of the terms in the stability condition [15] by inserting 
typical values of the parameters involved: For example, 
consider a 5000-lb thrust rocket airplane installation: 


Length of tubing from tanks to motor = 5 ft 


Equivalent diameter of tubing = 2.5 in. 

Mass flow of propellant = 25 lb/see = 0.78 slug/sec 
Chamber pressure = 400 psia 
Tank pressure = 550 psia 
Characteristic velocity (c*) = 5000 fps 


Gas constant R of chamber gas = 2000 ft-Ib/slug, R 
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Adiabatic flame temperature = 5500 R 
Characteristic length of motor (L*) = 60 in. 
_ 5 X 0.78 5000 X 5 300 


= 400 x 4.9 T 2000 x 5500 409 ~ °-002 + 0.002 = 0.004 see 


Ter 

In this case, the system will be stable if the combus- 
tion time lag is less than 4 milliseconds. On the other 
hand, if the same motor is installed in a similar system, 
but with 30 ft of tubing, the tolerable limit for the time 
lag is increased to 14 milliseconds. A propellant com- 
bination and combustion chamber that chugs in the 
first installation might be stable in the latter. 

It is not unreasonable to expect the actual time 
lags of practical propellant systems to fall within this 
range. Although the process of ignition is not quite 
like that of combustion, the observed ignition lag of 
self-igniting propellants may be expected to indicate 
an upper estimate of combustion time lags. Observed 
ignition lags, under conditions similar to the injection 
process in a rocket motor, are of the order of 5 to 30 
milliseconds, i.e., in the sensitive range (3). 

An interesting interpretation of condition [15] can 
be developed by considering the over-all combustion 
reaction. The time lag 7 may be defined as the interval 
between the instant of entry of an elementary mass of 
liquid propellant and the later ‘instant’? when this 
mass is converted to high-temperature gas, that is, 
when it exerts its full effect upon a pressure gage con- 
nected to the combustion chamber. This time interval 
7 includes the time required (on the average) for the 
gasification of the elementary liquid mass and the time 
required for the gas-phase reaction. The gasification 
time is not calculable on simple grounds, particularly 
since the process of vaporization is always accompanied 
by gas-phase or liquid-phase reactions of unknown 
character. The gas-phase reaction time, on the other 
hand, can be estimated (roughly) by equating it to 
the gas-phase residence time in a combustion chamber 
of minimum volume, i.e., an empirically determined 
volume such that further reduction results in a serious 
loss in specific impulse. This residence time is given 


by: 
( Ve)min Pa(mean) c*L* min 
/Po(mean) Po(mean) Po(min) RT 
A reasonable assumption is that 
then, 
min 
Ty RT. [ ] 


Denoting the time in the liquid phase as 7,, condition 
[15] takes the following form as a condition for stabil- 
ity: 

lin 2Ap 
PAs RT. RT, 

The second term on the left side and the correspond- 
ing term on the right side are usually of the same order 
of magnitude, since, in practical designs, 2Ap/p, is 
less than unity and L* is generally made greater than 
L* nin. There are sound reasons for such design 
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practice (based upon considerations other than the 
control of chugging); therefore, the inequality becomes 
sensitive to the relative magnitudes of the term on the 
left describing the momentum of the liquid in the feed 
lines and the term 7, on the right. 

This comparison indicates that in the case of self- 
igniting propellants (e.g., nitric acid and hydrazine) 
chugging is less likely to occur than in the case of non- 
self-igniting propellants (e.g., nitric acid—kerosene) by 
virtue of the relatively small gasification time of the 
former. In fact, chugging may be impossible with 
certain self-igniting combinations even if I. is very 
small, since the L* term may insure the inequality. 
On the other hand, care may be necessary in rockets 
burning less reactive liquids to provide sufficient 
momentum in the feed lines. 

Condition [15] indicates that, in the trend toward 
compact liquid rocket installations with short pipe 
lines, the appearance of instability will become more 
probable. Of the remedies available to the designer 
(listed above), those involving changes in the propel- 
lant feed system are simpler than those involving 
changes in the configuration of the rocket motor. The 
pressure difference Ap may be increased by inserting 
flow restrictors in the lines, but this has the disadvan- 
tage of requiring a higher supply pressure. A more 
attractive alternative is to increase the length or 
decrease the area of the feed lines: This would require 
a smaller increase in supply pressure to achieve the 
same degree of stability. 

In computing the magnitude of the momentum term, 
if the feed system consists of several passages in series, 


that term in condition [15] becomes 3 In such 
cs i 

calculations, the momentum of the liquid in the re- 
generative cooling passages should not be overlooked. 
In systems having centrifugal pumps instead of pres- 
surized propellant tanks, the theory is directly ap- 
plicable, but the momentum term should include the 
suction lines as well as the high-pressure lines. In the 
case of a system with positive displacement pumps, 
this theory would deny the possibility of unstable 
flow (unless a compressible air cavity or elastic tubing 
is present to permit flow oscillations). 

In systems with long feed lines it is possible that the 
effective length may be somewhat less than the meas- 
ured length due to compressibility of the liquid, par- 
ticularly if bubbles are present. A rough calculation 
indicates that the maximum effective length is of the 
order of (m/AzpBpw), where 6 is the volume com- 
pressibility of the liquid. At 100 cycles per sec, this 
would be of the order of 10 to 20 ft for the usual steel 
tubing filled with air-free nitric acid. (However, a 
more precise analysis including the role of the elasticity 
of the liquid should be carried out for such cases.) 

Caution is required in the application of these 
stability criteria to cases in which the linearization 
performed herein may be invalid. In addition to the 
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mathematical linearization above, a physical lineariza- 
tion is tacitly assumed in the statement that 7 is a 
constant independent of pressure. 

In conclusion, it is suggested that experiments be 
conducted to test the correctness of this theory; first, 
to determine the magnitude of the hypothetical time 
lag and its dependence on pressure, injector configura- 
tion, etc.; and second, to determine whether conditions 
[15] and [19] are valid conditions for stability. 
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Single Flow Jet Engines—A Generalized Treatment 


By J. V. FOA 
Head, Propulsion Branch, Cornell Aeronautical Laboratory, Inc , Buffalo, N. Y. 


Generalized performance equations are obtained for 
single flow air breathing jet engines from entropy considera- 
tions. The analysis is simplified by the use of convenient 
parameters in the description of the modes of compres- 
sion, combustion, and expansion. The generalized equa- 
tions are applied to a few cases of current interest, for the 
purpose of illustration. 


IR-BREATHING jet engines may be classified 
into two groups: (1) Single-flow machines, in 
which all the working medium undergoes the same 
heating process; and (2) multiple-flow machines, in 
which two or more portions of the working medium 
undergo partially or wholly different processes. 

The purpose of this paper is to present a generalized 
method of performance analysis of existing or conceiva- 
ble power plants of the first group. The method is 
based on the description of the modes of compression, 
heating, and expansion of the working fluid in terms 
of parameters that can be varied continuously over 
their respective ranges of practical significance. The 
thrust and specific impulse of single-flow jet engines 
are then derived, in the most general form, as functions 
of such parameters. 

This method yields a direct comparison of the per- 
formance of existing single-flow jet engines, together 
with an indication of the extent to which they may 
be improved and of what may be expected of other 
conceivable types of engines of the same group. 


Symbols 


= speed of sound, fps 

= cross-sectional area of duct, sq ft 

= specific heat at constant pressure, Btu per lb, ° R 

Cy = specific heat at constant volume, Btu per Ib, ° R 
= frequency of pulsations, cps 


ope 
I 


g = acceleration due to gravity, ft per sec? 

H = heating value of fuel, Btu per lb 

I = impulse, Ib sec 

I, = air specific impulse (impulse per Ib of air), sec 
I, = impulse per cycle, lb sec 

I; = I,/g = fuel specific impulse (impulse per Ib of fuel), sec 


= mechanical equivalent of heat 
= mass flow of air through the engine, slugs per sec 


M = Mach number 

n = exponent of the polytropic describing the heating 
mode 

p = static pressure, lb per sq ft absolute 

P = stagnation pressure, lb per sq ft absolute 

q = heat added per |b of air, Btu per Ib 

{a = total amount of heat added per lb of air, Btu per lb 
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= total amount of heat rejected per lb of air, Btu per 


dr 
Ib 
r = compression ‘atio through compressor 
R = gas constant, ft-lb per Ib, ° R 
8 = specific entropy — tu per lb 
t = time, sec 
= thrust, lb 
u = velocity, fps 
v = specific volume, cu ft per lb 
V = volume of heated region, cu ft 
Y = ratio of specific heats 
Ye) Yhy Ye = Specific heats under various temperature conditions 
Ne = polytropic compressor efficiency, defined through 


= 1 
the equation r = ? 


nh = combustion efficiency 
= polytropic turbine efficiency, defined through the 


ne 
"(ve — D 
equation P, 
v = static temperature, ° R 
0 = stagnation temperature, ° R 
r = product of all downstream-to-upstream stagnation 
pressure ratios associated with internal drag 
p = density, slugs per cu ft 
o = product of all downstream-to-upstream stagnation 
pressure ratios associated with shocks 
¢ = fuel/air ratio 
Subscripts 
0 = free stream 
1 = downstream end of inlet diffuser; compressor inlet 
2 = entrance to heating region 
3 = end of heat addition; turbine inlet 
4 = turbine outlet 
e = exhaust (at tailpipe exit) 


Assumptions and Approximations 


1 M;, and M, are assumed to be low enough that 6, ~ 
and 0, = ve. 

2 It is stipulated that the amount of heat added and the 
mechanical work performed by the gas during the heating process 
be always positive or zero (which means that the specific volume 
of the fluid will never decrease during the heating process). 

3 Heating is assumed to involve no change in the chemical 
composition of the working fluid. 

4 All particles of the working fluid are assumed to undergo 
the same thermal cycle in each engine considered (except for the 
purging of residual gases). 

5 Changes in the mass rate of flow of the fluid passing 
through the engine are neglected, and the working fluid is as- 
sumed to be air at all sections of the flow path. 

6 The effect of temperature on the ratio of specific heats for 
air is taken into account by the use of mean effective values, as 
follows: 


during compression y = ye = 7/5(¢pc = 0.24, Cve = 0.171) 
during combustion y = ya = 1.36(cpr = 0.259, con = 0.19) 
during expansion y = ye = 4/3(cpe = 0.275, Cue = 0.206) 
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ne : 
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7 Heat losses through the walls of the engine are neglected 
in all cases. 

8 The entire work output of the turbine is assumed to be 
utilized to drive the compressor. 

9 Complete expansion to ambient pressure within the exhaust 
nozzle is assumed in all cases. In the case of nonsteady-flow 
engines, where flow conditions change periodically in a nozzle of 
fixed configuration, this assumption has no meaning as such, 
except in the subsonic phase of the discharge, and must be inter- 
preted as an approximation. With this assumption, the thrust, 
defined as the change of momentum of the air which flows 
through the engine in the unit of time, is 7’) = m(ue — wo) in 
the case of steady-flow engines, and 7) = m(a, — wo), where 


f [m/f ; 
te = =I: udm, in the case of pulsating-flow engines. 


Representation of the Heating Mode 


Let the change of state which accompanies heat 
addition be approximated by the polytropic relation 
po” = constant. This representation is no doubt 
open to challenge, but so little is known about the 
actual process of heat addition by combustion that 
there seems to be little ground for preference in the 
matter of approximation, so long as the underlying 
assumptions appear to be tenable in the light of current 
experience. The advantage of the polytropic approxi- 
mation is that it provides a parameter—the polytropic 
exponent n—which may be varied continuously over 
a wide range and is alone sufficient to describe an in- 
finite variety of heating modes of practical interest. 

The First Law of Thermodynamics yields 


On the other hand, assumption 2 may be expressed in 
the form 


dq/d# > Coh 
which is compatible with Equation [1] only if 
n < 1.0. It follows that the range of significance of n 


in this analysis is from — © (isochorie mode) to +1.0 
(isothermal mode). 

From Equation [1] and the polytropic relation 
dp/d3 = np/(n — 1)8 one obtains 


dp_p _\>Owhenn<Oorn> (2) 


Equation [2] shows that, within the range of signifi- 
cance of n in this analysis, the static pressure of the 
gas decreases as heat is added to it when n is positive 
and increases when n is negative. Concurrently with 
this pressure variation, the flow velocity will increase 
as heat is added if 0 < n <.1.0 and decrease if n < 0. 
The extent to which the isochoric mode can be ap- 
proached without becoming incompatible with the 
maintenance of stationary flow in the heating region 
depends primarily on the shape of the heating duct 
and on the distribution of heat sources within it, but 
is very limited in any case. Indeed, Equation [1] 
yields 


(n — 1)qa 


(3) 


which, within the stipulated range of interest of n, is 
compatible with the stationary flow condition 
only if 
Av, — 1) qa 
(ye — 1) Me? 


This inequality expresses a necessary (but not 
sufficient) condition for the maintenance of stationary 
flow in the heating region during the heating process. 
The right-hand member of the inequality is the value 
of n for which the addition of a dimensionless quantity 
of heat ¢,/Cp,3%2, would cause the gas to come to rest 
from an initial Mach number Mp». It represents, 
therefore, the limiting value of » for stationary flow. 
It can be verified that this value is always very close 
to zero. Therefore, the line of demarcation between 
stationary and nonstationary processes can, for all 
practical purposes, be set at n = 0 (isobarie mode), 

Equation [2] can now be interpreted to mean that 
stationary flow can be maintained when n is positive 
and that heating processes in which n is negative will 
generally result in pressure fluctuations, the amplitude 
of which will increase with the absolute value of n, 
It has been shown (1)! that, in the nonstationary range, 
the absolute value of n is—all other conditions being 
equal—directly proportional to the rate of heat release. 
The heating mode may, however, be largely affected 
also by other causes, such as pressure waves or dis- 
sipative agents. One may, indeed, also regard as valid 
the converse interpretation of Equation [2]; -namely, 
that, if the static pressure increases during the heating 
process, then n is negative. In other words, just as the 
heating mode may be responsible for the generation 
of flow fluctuations, so may flow fluctuations which 
are otherwise excited be responsible for modifications 
of the heating mode. The parameter n is used here to 
describe the changes of state which accompany the 
addition of heat, regardless of whether they are pro- 
duced by heating alone or also by other causes. 

The present analysis will accordingly consider two dis- 
tinct classes of jet engines, corresponding to two 
distinct ranges of the parameter n: steady-flow engines 
which will be analyzed for values of n ranging from 
n = | (isothermal heating) to n = 0 (isobaric heating); 
and pulsating-flow engines, the analysis of which will 
cover all polytropic heating modes isobaric 
(n = OQ) to isochoric (n = —~), 


Exhaust Temperature 


From assumption 4, it follows that all the working 
fluid (with the possible exception of “residual’’ gases 
in pulsating-flow engines) undergoes the same total 


'‘ Numbers in parentheses refer to references on page 126. 
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specific entropy rise (s, — 89) in flowing through the 


engine. On the other hand, because of assumption 
9, 
— 8 = 
where 
(43/9) 
Assuming, for the sake of simplicity, ¢ = ¢,,, one 
obtains 
Se — 


Entropy Rise Through the Engine. Thermal 


Efficiency 


The specific entropy rise (s, So) will now be cal- 
culated as the sum of the following increments: 
(a) Entropy rise through shocks: 


R 
Ass = — —Ina...... 


where o is the product of all the downstream-to-up- 
streim stagnation pressure ratios across the various 
shocks which are traversed by a particle of the working 
fluid as it flows through the engine. 

(b) Entropy rise due to internal drag: The effect 
of wall friction and flow separation in subsonic dif- 
fusers, sharp changes of cross section, narrow passages 
and bends, or around such obstacles as fuel injectors, 
flame holders, ete., is in each case a loss of stagnation 
pressure without change of stagnation temperature. 
The entropy rise due to internal drag is therefore 


where A is the product of all the downstream-to- 
upstream stagnation pressure ratios associated with 
such losses within the engine. 

(c) Entropy rise through compressor: 


Ase = Cpe In 
0. 
= Cpe 1 — ne) In [7] 
1 
or 


(d) Entropy rise due to heat addition: Equation [1] 
yields 


whence 


Va . 
4a, © li = 
—1) 


Finally, from Equation [3], and with assumption 1, 
one obtains 
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= 
in = |S] 


It can be verified that, when As, is expressed in this 

form, = (As,) > 0: with a fixed value of g,/€pn42, As, 
n 


increases monotonically as n is increased from — © 
to +1.0. 

In dealing with engines in which material properties 
impose a limitation on the peak stagnation temperature 
of the cycle, it is convenient to express As, in terms of 
8s. 

In the stationary flow range (0 <n < 1), one has 


6. qa 
C [9] 


From Equations [8] and [9], one obtains 


n— n(n — 1) 


y(n — 1) a ( _ a ) 
(n ya) 


and it can be verified that, in this range, > (As,) > 0. 
n 


As; = Cpa 


In the nonstationary flow range (n < 0), Equation 
[9] no longer applies, but .W3 < Mo, and it is therefore 
permissible, in accordance with assumption 1, to sub- 
stitute 6; for 33 as well as @ for d in Equation [3]. 
This equation then becomes a relation between 63; and 
6. which, together with Equation [8], yields 
In 


(n — ya) 


_ 


As, = —C, 


and in this case one finds > (As,) < 0. In other 
n 


words, when g_/Cpr93 is fixed, ncn’) is negative in the 
n 


nonstationary flow range and positive in the stationary 

flow range. It follows that the minimum value of 

the entropy increment As, with the prescribed value 

Of qa/Cpx9s is obtained in the neighborhood of n = 0. 
(e) Entropy rise through the turbine: 


AS: = Cpe In 7 In P, 
n 65 


The entropy change due to chemical reaction is 
neglected in this analysis. The resulting error is of 
the same order of magnitude as that involved in the 
assumption of constant specific heats during the com- 
bustion process. At any rate, the error depends 
primarily on the fuel used and will not appreciably 
affect the relative merits of different types of engines 
employing the same fuel. 

As, is almost always the largest of all entropy incre- 
ments within the engine, and it is on the magnitude 
of this increment that the thermal efficiency of the 
engine primarily depends. In accordance with as- 
sumptions 7 and 9, the amount of heat rejected per 
unit weight of the working fluid in each cycle is g, = 
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€,(8 — &%). Therefore, from the definition of thermal 
efficiency and Equation [4], one obtains 


From inspection of Equation [11] and from the dis- 
cussion of Equations [8], [8’], and [8”], it may now 
be concluded that, when all entropy increments except 
As; are assumed to be constant, the thermal efficiency 
of a jet engine may be affected by changes of heating 
mode in two distinctly different manners, depending 
on whether OF Ga/Cpx9s is prescribed: in the 
former case, n, increases monotonically as the heating 
process approaches the isochoric mode; in the latter 
case it becomes maximum when the heating mode is 
approximately isobaric. 


Propulsive Efficiency 


With the assumption of complete expansion of the 
working fluid at the exhaust (p, = go) the thermal 
and propulsive efficiencies can be expressed in the 


forms 
pat 
J, ( 


fe (ue — 
(ue? — uo?)dm 


respectively. It is clear that the type of discharge 
which is most favorable from the standpoint of pro- 
pulsive efficiency, for a given thermal efficiency and 
for given uo, m/f, and qa, is that described by the func- 
jie = u.(m) which maximizes the flux 
udm through the exit, with the constraint u,2dm 
= constant. Euler’s equation for the vanishing of 
the first variation of the integral i 4 Ue — Au,?)dm, 
where A is an “arbitrary” constant to be determined 
from the equation of constraint, yields immediately 

= 0.54 = constant as the extremal, indicating 
that the maximum propulsive efficiency is obtained 
when all the working fluid leaves the engine at the 
same velocity. It must be noted that this extremal 
may be interpreted to represent either a stationary 
flow at the exhaust or an exhaust velocity which is a 
“square wave’ function of time—the gas being dis- 
charged all at the same velocity but intermittently 
with respect to time. 


Ue? — uUo?)dm 


and 


Np 


Stationary Flow Engines 


Assuming again, for the sake of simplicity, that the 
effective mean value of c, through the engine be é, = 
Cp, (a simplifying assumption which is certainly justified 
by the degree of approximation of this analysis), the 
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exhaust velocity of a stationary-flow engine can be 
calculated from the energy relation 


Ue? = Uc? + cpr (# 
which, together with Equation [4], yields 
Ue? = ug? + epr oof exp ( ) 


and therefore 


: if Ue — Uo 


= 
mg g 


Se — 80 


With the values of y and c, assumed in this analysis, 
one obtains 


}12.500] 1 + 0.185M.? + 


qa Se — So '/2 


If heating takes place by combustion, gq = H¢nn; 
hence, J; = 
qa 


Pulsating Flow Engines 


The generation of thrust in intermittent jet engines 
is a periodic or almost periodic process of considerable 
complexity, because each cycle is affected by the 
waves which still reverberate through the engine from 
the previous cycles and similarly affects the cycles 
which follow. The method of characteristics can be 
applied to the investigation of the approach to periodic 
operation, i.e., to such conditions as are finally estab- 
lished when the effect of previous history has become 
the same for all cycles. A large accumulation of 
errors is likely to occur, however, in the graphical 
construction of the extensive characteristic network 
which is required to arrive at a description of the 
periodic process. While it may be argued that the 
accuracy of the results is still within the approximation 
of the conjectural assumptions which must be made at 
any rate with regard to the combustion mode, it is 
clear that with the very use of such approximations 
the laboriousness of the method becomes unwarranted, 
regardless of its accuracy. Indeed, there is reason to 
believe that a comparable degree of over-all accuracy 
could be obtained by the use of assumptions and ap- 
proximations which will permit direct treatment of the 
periodic operation. It is, e.g., more convenient (and 
perhaps more accurate) to apply the polytropic ap- 
proximation to what, on the basis of observations, may 
reasonably be expected to be the combustion process in 
periodic operation than it is to apply it to the initial 
explosion only, with an additional set of assumptions 
as to the effect of pressure waves on the heating mode. 
It must be noted that the waves which propagate 
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through the heating region when each explosion takes 
place and those which are generated in each explosion 
are, in general, of predominant importance in the 
economy of -the- eyele.-.-The. simplification. which: is 
brought about by the use of the polytropic approxima- 
tion is particularly important with regard to these 
waves, because the parameter n describes the changes 
of state which occur during the combustion process 
regardless of their genesis. Therefore, all waves which 
propagate through the heating region are implicitly 
accounted for not only in so far as they modify the 
heating mode but also as they set out to affect the 
succeeding phases of the cycle, and previous history 
need not be given further consideration. 

The thrust of intermittent jet engines can then be 
assumed to be obtained from the periodic discharge 
of gases from regions where the pressure is periodically 
raised by explosions, the frequency of the forcing 
function being approximately the same as the damped 
natural frequency of the oscillating gas system on 
which it acts. 

With these assumptions, the changes of state or flow 
parameters at any station of the flow path during one 
cycle of the periodic operation become identified with 
the corresponding changes during the first complete 
fluctuation in the free oscillating discharge which 
would follow a single explosion of the type considered 
in the periodic process. In particular, the impulse 
which is generated in each cycle can be obtained from 
consideration of this single fluctuation, the initial 
condition being the same as that which is produced, 
in the periodic operation, by each explosion. 

The analysis by the method of characteristics is 
thus reduced to computations over a single cycle. The 
impulse generated in each cycle can then be calculated 
through time-integration of the momentum transport 
at the exhaust, for any given free stream velocity: 


I. = — uo)dm 
= J, PeActle( — Uo)dt 


where 7 is the time of outflow of the mass m/f. 

It must be noted that 7 does not represent the entire 
period 1/f of the exhaust velocity fluctuation but only 
that portion of it during which the exhaust velocity 
is positive. Indeed, if m/f represents all the mass of 
air which flows into the engine and is brought to rest 
there in each cycle, then the backflow impulse con- 
tribution is already accounted for in the momentum 


term — 7 and if, as is common practice, m/f is 
taken to represent only the mass of air which is heated 
in each cycle, then the backflow air is not part of this 
mass and must be given separate consideration. 

The approximations introduced so far simplify the 
problem to a considerable extent but are not sufficient 
to make it amenable to analytic solution. In particu- 
lar, there still remains the important question as to 
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what effect changes of shape of the reservoir and of the 
discharge duct may have on the impulse which is 
generated in the discharge of a given mass of gas from 
given initial reservoir conditions to a given exhaust 
pressure. What can be said by way of an answer is 
still tentative. A rather extensive investigation of 
blow-down discharges, which has recently been carried 


_ out at the Cornell Aeronautical Laboratory by the 


method of characteristics, has revealed that, under the 
conditions stipulated above, changes of shape have 
important effects on the duration and detailed history 
of the blow-down discharge, but hardly any discernible 
effect on the total impulse which is generated thereby. 
This observation suggests that the impulse per cycle, 
with any given reservoir volume, exhaust pressure, 
and initial reservoir conditions, may—at least in the 
absence of important entropy gradients in the discharge 
duct—be substantially independent of the shape of the 
duct and therefore also of the duration of the blow- 
down discharge. 

Let this conjecture be used as a working hypothesis. 

Now, if the discharge time is very long compared 
to the travel time of waves over the length of the duct, 
the exhaust velocity can be assumed to be at every 
instant the same as it would be in the stationary isen- 
tropic expansion from the same instantaneous reser- 
voir condition to the same exhaust pressure; in other 
words, wave phenomena can be neglected and the dis- 
charge can be treated as a quasi-stationary process. 
Therefore, according to our hypothesis, the impulse 
generated in the blow-down discharge of a gas from a 
container and through a duct of any shape should be 
approximately equal to that which would be generated 
in the quasi-stationary discharge of the same mass of 
the same gas from the same initial reservoir conditions 
to the same exhaust pressure. 

The operation of intermittent jet engines has been 
analyzed as a quasi-stationary process by several 
investigators (2,3). It is interesting to test the validity 
of such approximations—and therefore also of our 
working hypothesis—by comparison of the results of 
the quasi-stationary discharge analysis with those 
obtained by the method of characteristics. 

As pointed out above, the expansion ends, as far as 
the generation of impulse is concerned, with the first 
vanishing of the exhaust velocity. In the analysis 
of the discharge as a quasi-stationary process, this 
condition is found to occur when the pressure in the 
chamber becomes equal to the ambient static pressure 
Po, i.e., in the absence of entropy gradients in the 
discharge flow, when the temperature and density 
throughout the chamber become equal to the constant 
exhaust temperature #, and density p,, respectively. 
Then, letting subscripts 3 and c denote initial and 
instantaneous reservoir conditions, respectively, and 
V the volume occupied by the heated gases when the 
expansion starts, the impulse per cycle is 
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I, = Ue dm dd. = V(ps Pe )Uo 


the residual mass Vp, being scavenged at free-stream 
velocity,? so that no impulse contribution is derived 
from its discharge. Now, 

ve 


— lds, 


P3 
—1 1) 


dm = —V-dp. = — 


and 


Ye 


Therefore 
Vos 


With the assumed value 4/3 for y,, one obtains 
I. = 0.1368Vpsa [s (5) + x 
3 
The value of the impulse calculated by means of 
this equation agrees with the above-mentioned results 
obtained by the method of characteristics for ducts of 
various shapes, within the accuracy of these results. 
These tests of the validity of our hypothesis are, of 
course, far from conclusive. While there seems to be 
little doubt that the quasi-stationary approximation 
can be used with confidence in the evaluation of certain 
items of performance, its validity should be supported, 


if possible, by more rigorous argument. 
The air specific impulse is then 


I, = 


10.1368 [s 129 15 | x 


again with g, = Hgm and therefore I, = I,Hm,/qa, if 
heating takes place by combustion. 
The ratio 33/3 is obtained from Equation [3]: 


Bo De Oy Do Con — 80 
and the ratio 3,/83 is calculated from Equations [4] to 
[8], [10], and [15], which, in the case of compressorless 
pulsating-flow engines (As, = As, = 0), yield 


| — 1) 
— n) 


=(or) 


_(n — 
+ — n)d2 


The beneficial effect of an increase of the ratio 


. [16] 


2 Failure to account for the cost of scavenging led to er- 
roneous results in reference (1). The author is indebted to R. 
Weatherston for calling his attention to this omission. 
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32/0, i.e., of precompression of the charge beyond 
the ram level, can be evaluated by means of these 
equations. It must be noted that it is only in non- 
stationary flow that such precompression can be pro- 
duced without the aid of a mechanical compressor, 
This possibility must indeed be regarded as one of the 
salient features of pulsating-flow jet engines. 

The “hammer effect,’ resulting from rapid inter- 
ruption of a flow, illustrates this possibility. Let the 
initial conditions of a flow in a straight tube be u = 
u; and a = a, Interruption of the flow by rapid 
closure of a gate generates an upstream propagating 
pressure wave (hammer wave). Let subscript H 
refer to the flow region between this wave and the gute. 
Then, uz = 0. The compression ratio produced in 
this manner, assuming isentropic changes of stute, 
can be computed from the relation (4) 


ay — a; = (us — un) = 2 


whence 


PH _ (1 re 
Pi 2 


4 
It can be verified that py > P,; if M; < “eae The 


27 
Mi)? 


hammer compression effect is of course lower where 
the hammer wave has steepened to form a shock. 

Another way of producing precompression in pul- 
sating-flow engines is by utilization of the waves which 
are generated in the combustion process. Equation 
[2] reveals that a nonstationary polytropic heating 
process is always accompanied by a pressure rise; and 
that this rise increases monotonically as the heating 
process approaches the isochoric mode. Therefore, 
intermittent combustion is a source of pressure waves 
which can be used to precompress the fresh charge. 

The manner in which pressure waves may be timed 
and utilized to precompress a charge depends in the 
first place on the mode of combustion. In the pre- 
ceding paragraphs, combustion was assumed to take 
place uniformly throughout the burning region—a 
condition which is perhaps approached in actual opera- 
tion when ignition takes place more or less simul- 
taneously at several points in the combustion chamber, 
or when the flame front is broken up by turbulence 
into a number of scattered burning “islands.” In 
this case, the pressure is found to rise uniformly in 
the burning gas, generating pressure waves which are 
propagated and reflected in the surrounding regions 
and may be utilized to precompress the fresh charge in 
the following cycle. 

On the other hand, any departure from the condi- 
tion of uniformity assumed in the preceding paragraphs 
will introduce a certain amount of precompression 
within the cycle, because the pressure waves generated 
by the burning portions of the mixture will precompress 
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the portions that are yet to be burned. An extreme 
case of wave precompression of the latter type is that 
which occurs in the presence of a plane normal flame 
front propagating axially through a duct filled with 
a combustible mixture. In this case, and assuming 
that the change of state across the flame be isobaric, 
the wave precompression ratio is found to be (1) 


where U is the ratio between the rate of flame propaga- 
tion relative to the unburned gas and the velocity of 
sound in the unburned gas. Utilization of this mode 
of precompression was recently attempted by Kahane 
in the “intermittent ramjet’’ (5). 

Equations [14] to [16] can be used to calculate the 
air and fuel specific impulse of a large class of single 
flow intermittent-jet engines. There exists, however, 
a class of engines in which the thrust-generating dis- 
charge departs considerably from the analytical model 
used in the derivation of these equations. This special 
class is that of the engines in which nonstationary 
processes are confined to those regions in which they 
can be utilized to increase the thermal efficiency, 
while substantially stationary-flow conditions are main- 
tained at the intake and exhaust, to reduce the intake 
shock losses and to obtain maximum _ propulsive 
efficiency (see page 118 under ‘“‘Propulsive Efficiency’’). 
The performance of any engine of this type must again 
be caleulated by Equation [12], which applies to any 
engine in which both inflow and outflow are stationary, 
regardless of the character of the flow within the 
engine. 


Examples 


The generalized theory is applied in this section for 
the purpose of illustration to the following families 
of power plants: (a) ramjets and (b) turbojets, with 
any stationary-flow heating mode; (c) pulsejets, with 
any nonstationary-flow heating mode; and (d) wave 
engines, with isobaric combustion. 

The results of the computations are summarized by 
charts for each of the engines considered, with the 
fuel specific impulse as the ordinate and the air specific 
impulse as the abscissa. Lines of constant air/fuel 
ratio (1/¢g) are also shown in the charts. These are 
straight lines passing through the origin, because 
I; =I a/ 

The following sets of numerical values, in addition 
to those listed under assumption 6, were used in all 
performance calculations: } = 0.9, no = 0.85, m = 
0.80, m = 0.90, H = 19,000 Btu/lb, and & = 416° R 
(corresponding to a9 = 1000 fps and to an altitude of 
about 29,000 ft in NACA standard atmosphere). 

The transition from supersonic to subsonic flow in 
the inlet diffuser was assumed to occur through a single 
hormal shock, and in all cases the diffuser configuration 
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FIG. 1. FUEL SPECIFIC IMPULSE AND AIR SPECIFIC IMPULSE OF 


RAMJET 


Was assumed to be that described in (6). In station- 
ary-flow engines the shock was assumed to be located 
at the throat of the diffuser. In nonstationary-flow 
engines the moving shock was assumed to be retained 
within the diverging portion of the diffuser, and its 
mean effective strength was assumed to be that of a 
normal shock occurring at the free-stream Mach 
number. 

The Ramjet: The ramjet can be broadly defined as 
any stationary flow (n > 0) compressorless (As, = 
As, = 0; 6. = 6) air-breathing jet engine. Accord- 
ingly, 

— 1 n= VA 


The performance of the ramjet can be calculated by 
means of Equations [12] and [18] for any set of values 
of n, qa, 8, and Mo. The results of such calculations 
are presented in Fig. 1. The performance of the ram- 
jet is seen to depend te a considerable extent on the 
value of n, all other parameters being equal. 

It is interesting to determine where the ramjet with 
the commonly used constant-area burner would fit 
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in this chart. Stationary heating of a flow in a con- 
stant-area duct is not a polytropic process, but for the 
infinitesimal change of state at any point the function 
p(p) may be replaced by a polytropic tangent, in the 
form dp/dp = np/p, where n will of course vary con- 
tinuously along the heating region. Now, in stationary 
flow in a constant-area duct one has dp/p = —u-du 
and du/u = —dp/p, hence dp/dp = u?. Therefore, 
uv = np/p = na*/y, whence n = yM?. The Mach 
number in the combustion chamber increases from a 
very low upstream value M, to a downstream value 
M; which may be as high as 1.0. Concurrently, n 
will increase through the heating region from an up- 
stream value very close to zero, through the isothermal 
value 1.0 at M = 1/V y, to the isentropic value y 
at M = 1.0. It can be verified, however (7), that most 
of the entropy rise takes place in the low-Mach- 
number low-temperature region of the flow. There- 
fore, as far as the total entropy rise is concerned, a 
stationary heating process in a constant-area duct with 
alow inlet Mach number is equivalent to a polytropic 
process which departs only slightly from the isobaric 
(n = 0). 

Fig. 1 shows that the optimum heating mode is the 
isobaric at all flight Mach numbers. It is evident, 
then, that little room is left for improvement of ramjet 
performance through modifications of the-heating mode. 
In fact, a more nearly isobaric process would require 
a diverging combustion chamber, and the consequent 
increase of external drag could easily offset the slight 
thrust advantage. On the other hand, some improve- 
ment of ramjet performance can be obtained by 
reducing shock losses and internal drag and by increas- 
ing the combustion efficiency. Considering the pos- 
sibility of such improvements, it seems reasonable to 
regard the curves for n = 0 as fairly representative of 
the potential performance of ramjets. 

In general, the trend of the curves of potential ram- 
jet performance is similar to that of the performance 
of conventional ramjets as calculated by the usual 
methods of analysis. With extremely lean mixtures, 
the fixed losses (shock and friction) are predominant, 
therefore the fuel specific impulse increases as the energy 
input is increased. With richer mixtures, the entropy 
ie due to heat addition is predominant, and therefore 
in increase of heat input causes the fuel specific im- 
pulse to decrease, although the air specific impulse 
ontinues to increase. The maximum of J; occurs 
it very low values of J,, but the latter may be increased 
‘nsiderably without too much sacrifice in fuel econ- 
imy. The highest values of J; are obtained at flight 
Mach numbers of about 2, with the inlet diffuser 
‘nsidered in this analysis. The use of improved 
liffusers would have little effect on the performance 
it lower Mach numbers, but would increase somewhat 
te flight Mach number for maximum fuel economy. 
The Turbojet: In the turbojet, n > 0 and 6 = 4% 
‘in the ramjet, but 6 > 6. Assumption 8 estab- 


lishes the condition 
Cp O3 — 04) = Cpc(O2 — 


Therefore, 
1— Cpe 
exp *) = (or) ™ (A= 
Cph 
= vh 
[2 + (n — 1)qa 
— — (qa/Cpa)] 


. [19] 


cpe (1 — nt) 
_ Cpe [03 — (qa/epn) — Oo) nt 
Cpe 63 


The performance of the turbojet can be calculated by 
means of Equations [12] and [19]. Fig. 2 is a plot 
of the results of such calculations, in which the stag- 
nation temperature at the turbine inlet (63) was as- 
sumed to be 2500 R. This value is high in comparison 
with current practice but may well be attained and 
even exceeded in the near future, as a result of im- 
proved blade-cooling techniques. Then 


h 


and 


Neye 
B 2500 66,000 7° 


Therefore, in the (J,, J;) plane for any given flight 
Mach number, the loci of points corresponding to the 
same compression ratio are straight lines passing 
through the origin. The values of r are indicated in 
Fig. 2 alongside those of the air/fuel ratio on the 
constant —¢ lines. 

The lowest significant value of r is of course 1.0, 
since, when r = 1.0, the turbojet becomes a ramjet. 
On the high side, the establishment of any limitation 
on the value of r would be outside the scope of the 
present analysis. Therefore, no such limitation has 
been set in this calculation, and it will be noted that, 
as a result of this omission, the performance curves 
have been extended to regions where the specified 
conditions would impose compression ratios far above 
the values which are currently used. 

The best performance is obtained at all flight Mach 
numbers with n = 0, but the effect of small departures 
from the isobaric heating mode is relatively unimpor- 
tant and vanishes completely when the air/fuel ratio 
and the compression ratio are very high. Therefore, 
as in the case of the ramjet, little if any performance 
improvement could be derived in the turbojet from 
modifications of the heating mode. Also, the air/fuel 
ratio which is currently used in turbojets (about 60) 
appears to be a good compromise between those re- 
quired for maximum fuel economy (70 to 80) and for 
maximum specific thrust (about 40). The perform- 
ance superiority of the model considered in the cal- 
culations over present-day turbojets can therefore 
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2500 — 66,000 7 a 
2.975 
a 
: 
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be attributed almost entirely to the higher permissible 
temperature at the turbine inlet. The higher tem- 
perature is established through the use of higher com- 
pression ratios, without substantial change in mixture 
composition. The use of richer mixtures is advisable 
only when a high specific thrust is desired, even at a 
sacrifice of fuel economy, as in take-off or during 
acceleration in flight. Thus, e.g., the same compres- 
sion ratio (about 9) would be required for maximum 
specific thrust at My = 0 (take-off) and for minimum 
fuel consumption at My ~ 2; but the air/fuel ratio 
would have to be increased, during the acceleration, 
from about 40 at take-off to about 60 at the cruising 
Mach number. 

The Pulsejet: The pulsejet has been defined (8) as 
a compressorless pulsating-flow engine which “derives 
its thrust from a natural oscillation of a valved aero- 
duct.” The valveless pulsejet may accrodingly be 
described as deriving its thrust from a natural oscilla- 
tion of a valveless aeroduct. The principle of opera- 
tion is essentially the same for both types of power 
plant and their distinguishing feature, in the com- 
pressorless pulsating-flow class, in a complete or almost 
complete absence of wave precompression. These 
power plants will, therefore, be characterized in the 
present analysis by the condition 3 = 3, = %, which 
must be introduced in Equations [14] to [16] for the 
computation of their performance. 

It must be noted that As, is not only the entropy rise 
through the inlet diffuser shock, but also that associ- 
ated with the pressure waves which are generated by 
each explosion. Only shock waves are significant in 
this respect, and it should be borne in mind that this 
entropy increment must be considered only in so far as 
it results from wave phenomena occurring outside of 
the burning region, because the effect of pressure waves 
within this region is already accounted for in the 
description of the heating mode. Therefore, considera- 
tion must be given to the additional entropy term only 
if the heating mode is such as to produce pressure 
waves that are steep enough to form early shocks and, 
thus, to affect a significant portion of the nonburning 
region. 

In general, the entropy rise caused by these shocks 
is found to be insignificant unless the rate of heat 
release is very high, i.e., unless the absolute value of 
n is very large. In the limiting case of constant- 
volume combustion, n = —~, it has been shown (1) 
that this entropy increment is approximately equiva- 
lent to that occurring across a single normal shock 
propagating at a Mach number M = 3#;/9,. 

The calculated performance of the pulsejet is sum- 
marized in Fig. 3. The effect of changes of heating 
mode is pronounced only in the lower range of |n]. 

Since it has been shown (1) that |n| is almost linearly 
proportional to the rate of heat release, it follows that 
it should become decreasingly profitable to seek per- 
formance improvements by increasing the rate of heat 
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release. As a matter of fact, only at the lower Mach 
numbers is the most economical heating process very 
nearly isochoric: as the flight Mach number increases, 
the importance of the shock losses associated with 
high heating rates also increases and therefore the 
optimum |n| decreases. 

The chart for Mp = 0 is of limited significance 
because in static operation and at low flight speeds a 
considerable amount of cold air is drawn into the tail- 
pipe of the pulsejet in the backflow which follows the 
blow-down discharge. The propulsive efficiency is 
thereby increased because of the increased density of 
the tailpipe gases and of the consequent redistribution 
of the outflowing kinetic energy over a larger mass of 
gas. 

A gross evaluation of this increase of propulsive 
efficiency, which is called “internal augmentation,” 
can be made as follows. Let average values be indi- 
cated by bars. Then, if the energy transfer could be 
effected without losses, 7,2 would be inversely pro- 
portional to p, and therefore p,ti, would be approxi- 
mately directly proportional to V pp. The magnitude 
of this effect is, therefore, very appreciable. As the 
flight seeed increases the augmentation effect decreases 
rapidly, not only because backflow is reduced but also 
because the momentum lost by the augmentation air 
as it is brought to rest into the tailpipe from the sur- 
rounding stream is in any case proportional to the 
velocity of this stream. 

The chart for My = 0 in Fig. 3 must therefore be 
regarded as quite pessimistic. Indeed, the measured 
fuel specific impulse of some pulsejets in static opera- 
tion (9) has already exceeded the best theoretical 
values indicated by this chart. 

Assuming that the rate of heat release be in each 
case that corresponding to the optimum value of 2, 
comparison of Figs. 1 and 3 reveals the interesting 
fact that the potential superiority of the pulsejet over 
the ramjet is not limited to subsonic flight conditions 
but extends well into the supersonic domain for most 
practical values of J,. At the highest flight Mach 
numbers considered in this calculation, the performance 
of the pulsejet is still comparable, if not superior to 
that of the ramjet when the required air specific impulse 
is very large. It is also interesting to note (Fig. 3) 
that the rate of heat release required for optimum 
pulsejet performance at the highest Mach number is 
well within the range of possible improvement. This 
analysis is, of course, based on the assumption that the 
great difficulties which are now encountered in the 
high-speed operation of pulsejets can be surmounted. 
These difficulties are largely due to the inability of 
pulsejets to sustain ram pressure in the combustiol 
chamber during the charging phase of the cycle. Ther 
is reason to believe, however, that this situation cal 
be remedied, either by the introduction of a constrit- 
tion in the tailpipe (to “choke” the outflow over 4 
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longer portion of the cycle) or by shrouding the tail 
end to produce a local rise of static pressure; and per- 
haps also by other means. 

In Fig. 3 the static performance curve for n = 
—5 at sea level is also plotted in the chart for My) = 0, 
to illustrate the effect of altitude on specific impulse. 
This effect is seen to be quite small. 

Wave Engines: This designation is now customarily 
applied only to compressorless jet engines in which the 
fresh charge is precompressed by wave processes which 
are generated and/or controlled by timed valving (10). 

The entropy increment As, rates particular attention 
in this case because the process of wave precompression 
may involve shock waves of considerable strength. 
The losses associated with these shocks represent the 
cost of energy transfer from the hot gases to the fresh 
charge, the wave-engine counterpart of the turbine 
and compressor losses in turbojets. It must be re- 
membered that a train of weak pressure waves can be 
far more efficient as a compressor than a single strong 
shock producing the same pressure rise. Therefore, 
important improvements of precompression efficiency 
may be expected to result from refined techniques of 
shock-formation control. 

Now let the following conditions be assumed: (a) 
Wave precompression takes place through v shocks of 
equal strength; and (b) heat is added at constant 
pressure following wave precompression. Let 6; and 
o, denote the downstream-to-upstream static tem- 
perature ratio and stagnation pressure ratio, respec- 
tively, through each of the vy shocks. Then 


de 


so that 


In (82/0:) 
R In E In ] 
5 0 i 


In of 
where go is the stagnation pressure ratio across the 
stationary shock in the inlet diffuser; and 3 = 3; — 
Cyn 
Within the approximation of this analysis, 3; = 6, 
3. = 6, and 3; = 63. Therefore, 


h— 1 YR — 1 
h 


- Indi 190] 


Usually, the exhaust velocity from each element of a 
wave engine is substantially a square-wave function 
of time, and the engine is made up of a multiplicity 
of elements operating out of phase, so that the exhaust 
from the engine is substantially stationary. Similarly, 
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FIG. 4 FUEL SPECIFIC IMPULSE AND AIR SPECIFIC IMPULSE OF 
WAVE ENGINE 


the air intake in the various elements takes place in 
such close succession from the common inlet diffuser 
that the flow at the engine intake may be regarded as 
stationary. Therefore the performance equation for 
the wave engine is Equation [12], in which the entropy 
term is now obtained from Equation [20]. The cal- 
culated performance, for 6; = 4500° R and 6, = 1.5 
(which corresponds (11) to o; = 0.83), is plotted in 
Fig. 4. The potential performance of the wave engine 
is seen to be very promising, at least from the stand- 
point of fuel economy. 
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Approximate Calculations of Specific Heats for 
Polyatomic Gases 


By R. V. MEGHREBLIAN 


Guggenheim Fellow in Jet Propulsion, Guggenheim Jet Propulsion Center, California Institute of Technology, 
Pasadena, Calif. 


This paper presents a method for the approximate cal- 
culation of the ideal specific heats of polyatomic molecules 
occurring in rocket combustion processes. The method 
is based on the assignment of characteristic frequencies 
to individual interatomic bonds, the simplification con- 
sisting in the use of the same characteristic frequency for 
a particular bond irrespective of the over-all structure of 
the molecule in which it occurs. The method can be ap- 
plied to molecules for which precise spectroscopic data are 
not available. 


Introduction 

HE thermodynamic functions of gases can be 

calculated from spectroscopic data by using stand- 
ard statistical methods (1, 2).'. At moderate tempera- 
tures the existence of excited electronic states for poly- 
atomic molecules may be ignored. Furthermore, in- 
teractions between vibration and rotation make only 
small contributions to the total energy. Finally, the 
rotational motion of polyatomic gaseous molecules may 
be assumed to be fully excited for temperatures of in- 
terest in connection with calculations on rocket en- 
gines. 

It is well known that the specific heat at constant 
pressure C, of polyatomic gases, subject to the approxi- 
mations listed above, is given by either of the two rela- 
tions 

Received May 22, 1951. 

' Numbers in parentheses refer to references at end of paper. 


TABLE 1, 
Stretching Bending 
vibration vibration 
Bond (em~!) 
C—C (aliphatic) 910 650 
C—C (aromatic) 1500 600 
C=C 1200 910 
‘=C 2080 375 
C—H 3000 1050 
C—O 1030 1120 
C=0 1740 780 
C—N 1000 450 
C=N 2220 240 
C—Cl 650 260 
C—F 1050 1200 
C—Br 610 950 
C—I 530 880 


3n-5 
Cy =2R+R ... (linear molecules)... (1] 
2 
or 
3n-6 
C, =4R+R E(u;) ....(nonlinear molecules). .. . [2] 
j=1 


Here R represents the molar gas constant, n is the 
number of atoms in the molecule under discussion, and 
E(u;) represents the Einstein function (3) for the jth 
vibration frequency. The characteristic numbers wu; 
are determined from the normal vibration frequencies 
v, according to the relation 


uj = hv; /kT 


where h and k represent Planck’s and Boltzmann’s 
constants, respectively, and 7 is the absolute tempera- 
ture. 


Characteristic Vibration Frequencies and the 
Caiculation of Heat Capacities 


It is possible to obtain approximate values for the 
specific heats by utilizing the fact that characteristic 
vibration frequencies between two atoms are not usually 
strongly dependent on the atomic groups to which the 
vibrating atoms are attached (4). Thus, one can as- 
sociate with a given chemical bond two basically dif- 
ferent types of motion, viz., stretching and bending 


CHARACTERISTIC WAVE NUMBERS OF CHEMICAL BONDS 


Stretching Bending 
vibration vibration 
Bond (em~) (em~) 

H—O 3500 1350 
H—S 2570 860 
S=O 1250 520 
S—C 690 280 
N—N bond in NO 2220 590 
N—N bond in N;0. 280 500 
N—N 1000 900 
N—H 3300 1200 
N—O 1270 660 
N= 1470 650 
Si—H 2180 930 
si—O 1050 400 
B—H 2550 1150 
B—F 890 690 
B—N?® 


“ The information available on compounds with B—N bonds was inadequate for assigning average normal vibration frequencies. How- 
‘ver, calculations of Cp for borazole (B;N;He), based on the assignments 1350 cm~! and 500 cm~', have yielded reasonable agreements 


vith the more precise Cp data obtained by Crawford and Edsall (4). 
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vibrations. This correlation is useful (5) in spite of the 
fact that the normal vibration frequencies represent 
motions of the entire molecular structure rather than of 
two isolated atoms. 

Previous investigators (5) have compiled a list of 
characteristic frequencies for some of the important 
chemical bonds. However, these lists are inadequate 
for performing extensive calculations on propellants. 
Consequently an investigation was undertaken to ex- 
tend this information to additional bonds, and to revise 
the frequency assignments by using the best available 
spectroscopic data. A large number of different com- 
pounds (4) were examined in order to obtain reasonable 
values for the characteristic frequencies. 

The suggested assignments, in wave numbers, are 
listed in Table 1. The small and large values corre- 
spond, respectively, to the bending and stretching vibra- 
tions. The numerical values given in Table 1 do not 
differ greatly from the assignments of other investiga- 
tors for the chemical bonds for which previous assign- 
ments are available (5). 

Specific heats of polyatomic molecules may be cal- 
culated from Equations [1] or [2] by using the data of 
Table 1 together with an averaging procedure intro- 
duced by Bennewitz and Rossner (5), which has not 
been applied extensively by other investigators. . It is 
found, in general, that the sum of the stretching and 
bending frequencies will be fewer than the 3n - 6, or 
3n-5, normal vibration frequencies. This discrepancy 
is, of course, the result of the existence of complex 
group vibrations in the molecule. The suggested pro- 


TABLE 2. SAMPLE CALCULATIONS OF SPECIFIC HEATS OF GASEOUS 
COMPOUNDS 


Tempera- Cp Cp Cy 
ture (theo- (approxi- (experi- 
Molecule (C) retical) mate) mental) 
CH;0H 137 Sas 12.9 11.8 
N20, 27 15.86 15.63 ae, 
NH; 0 8.4 
300 ae 11.06 10.27 
500 pre 12.51 11.31 
700 ae 13.74 12.92 
850 Eve 14.52 13.22 
1000 15.15 14.65 
20 19.1 20 21.8” 
100 24.5 25.7 25.8 
350 38.35 39.45 38.9 
CH;NH, 10 11.72 11.85 
25 cae 12.29 12.41 
50 12.62 13.04 13.79 
C:He 64.5 13.6 14.1 13.44 
N20 20 10.24 10.14 9.24° 
(CH;),NH 10 15.8 15.8/ 
50 ity 18 18.66 
CHCI,Br 27 16.2 16.19 
427 21.48 21.72 
C.H;CN 25 15.24 15.04 
327 23.41 23.7 
* Cf. Ref. 6, p. 1228. 4 Cf. Ref. 7, p. 2316. 


” Cf. Ref. 7, p. 1275. 
© Cf. Ref. 7, p. 2317. 


cedure for evaluating C, is to account for the remaining 
unknown frequencies by assigning to the specific heat 
contribution of the undete:mined frequencies an aver- 
age value based on the contribution of the known bend- 
ing vibrations.2 Thus, the necessary adjustment is 
made by applying a correction factor (the ratio of the 
total number of nonstretching vibrations to the number 
of bending vibrations) to the C, contribution of the 
bending vibrations. Illustrative examples, showing 
the method of calculation, may be found in the p:iper 
of Bennewitz and Rossner (5). 

The results of sample calculations of specific heats 
are summarized in Table 2. Both theoretical (base: on 
actual frequency assignments) as well as experimental 
values of C, have been included. It is expected that 
results comparable to those indicated in Table 2 will be 
obtained for most polyatomic gases. For practical 
purposes, the calculated results should have acceptable 
accuracy. The principal advantage of the present 
method of calculation is the relative ease with which one 
can evaluate C, for polyatomic gases. 
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The Hypergolic Reaction of Dicyclopentadiene 
with White Fuming Nitric Acid 


By C. H. TRENT? and M. J. ZUCROW* 


The photographic study of the hypergolic reaction be- 
tween dicyclopentadiene and anhydrous nitric acid in- 
dicates that a solid phase is formed prior to ignition and 
that the ignition of the propellant is propagated from the 
solid phase. The latter has been isolated. Although the 
identity of the solid phase remains unknown, characteriza- 
tion tests have indicated that it is a nitrated polymer con- 
taining the nitro group and the nitrate ester group. 


HE discovery that dicyclopentadiene was hyper- 

volic with nitric acid led to the further study of 
that hydrocarbon for the purpose of obtaining informa- 
tion on the behavior of hypergolic hydrocarbon com- 
pounds. 


Photographic Study 


When liquid dicyclopentadiene and white fuming 
nitric acid were reacted in an inert atmosphere of dry 
nitrogen, spontaneous ignition occurred in spite of the 
absence of atmospheric oxygen. Furthermore, di- 
eyelopentadiene reacts hypergolically even when the 
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0.0035 sec after ignition 0.0056 see after ignition 
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PHOTOGRAPHS OF THE HYPERGOLIC REACTION BETWEEN DICYCLOPENTADIENE AND NITRIC ACID. 
FILM SPEED 1500 FRAMES PER SEC. 


nitric acid is added to it dropwise (approximately 10 
drops per sec) in an apparatus arranged so that a rapid 
stream of dry nitrogen sweeps the surface of the re- 
actants. If the nitric acid is added very slowly, there 
is formed, in the bottom of the reactor, a reddish-black 
solid substance which upon the further, more rapid 
addition of nitric acid bursts into flame. From these 
observations and others tobe discussed later, it appeared 
that a solid phase was formed prior to the ignition of 
the dicyclopentadiene-WFNA reaction. To establish 
whether or not a solid phase is formed, the reaction was 
studied photographically; a Western Electric Fastax 
high-speed motion-picture camera recorded the progress 
of the reaction. The photographs comprising Fig. 1 
are from portions of the film taken at 1500 frames per 
sec; that is, approximately 0.0007 sec between frames. 
The hypergolic reaction was produced by allowing 
fine streams of the reactants to impinge on a v-shaped 
metallic surface. To eliminate any effect of the presence 
of air on the reaction and to prevent the clouding of 
the metallic surface by the considerable quantity of 
vapor evolved during the reaction, the experiment de- 
picted in Fig. 1 was carried out under a rapidly moving 
atmosphere of dry nitrogen which swept across the 
surface of the reactants. The sequence in Fig. 1 clearly 
shows the formation of the solid phase (the dark area 
to the right of the injection apparatus) and that igni- 
tion occurs on the surface of that solid material. The 
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point of ignition can be seen as a faint light spot in the 
center of the solid (Fig. 1, frame 3) which grows in in- 
tensity in succeeding frames. Approximately 0.0056 sec 
after ignition has occurred, the combustion of the resid- 
ual vapor above the solid phase was initiated (Fig. 1, 
frame 7). A second point of ignition on the surface of 
the solid phase is seen to occur in Fig. 1, frame 8. 

Another series of photographs taken of the reaction 
between dicyclopentadiene and nitric acid resemble 
those of Fig. 1, differing only that the ignition occurs at 
a different place on the surface of the solid phase. If 
the reaction is carried out in still air, the ignition still 
occurs on the surface of the solid but the flame is rapidly 
propagated throughout the vapor above the surface of 
the solid. 

That the dark area on the metallic surface shown in 
Fig. 1 is actually a solid phase has been demonstrated 
by the actual isolation of solid compounds during con- 
trolled reactions between dicyclopentadiene and nitric 
acid. 


Chemical Nature of the Reaction 


The dropwise addition of a molar equivalent of nitric 
acid (99 per cent) to dicyclopentadiene immediately and 
exothermically forms an amorphous dark-red solid in- 
soluble in the original hydrocarbon. The solid is 
sufficiently stable to permit isolation and chemical 
examination. Tests showed that the solid compound 
or mixture of compounds formed are not shock sensitive, 
but when rapidly heated the solids burst into flame. 
Moreover, the addition of about 10 drops of concen- 
trated nitric acid to a small amount of the solids heated 
to 100 C will cause combustion to occur. 

The reaction between dicyclopentadiene and nitric 
acid has also been conducted in dilute solutions employ- 
ing carbon tetrachloride and petroleum ether as the 
diluents. From those reactions sufficient amounts of 
solid material were isolated for further study. It was 
found that the quantity of solids precipitated during 
the reaction increases as the molar ratio of nitric acid to 
dicyclopentadiene increases. The physical properties 
of the solid phase also appear to be altered, Table 1. 
In addition to solid compounds, very unstable liquid 
products were also isolated. The latter explode upon 
attempted distillation at atmospheric pressure. 


TABLE 1 SOLID FORMATION DURING REACTION BETWEEN 
DICYCLOPENTADIENE AND NITRIC ACID 


No. Dep,g HNO;,g Solidformed,g  D.p., C 


203 13.2 6.3 3.8 (I) 155-160 
204 13.2 "12.6 9.6 (VI) 172 
206 13.2 18.9 10.4 (XI) 210 


Insufficient research was conducted on the chemical 
properties of the solid compounds of Table 1 to estab- 
lish their exact chemical structure. The chemical be- 
havior of solids I and VI, however, leads to the con- 
clusion that the solid compounds contain the nitro 
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group and the nitrate ester group. The elemental 
analyses of I and VI are presented in Table 2. 


TABLE 2 ELEMENTAL ANALYSES OF SOLIDS DERIVED 
FROM REACTION BETWEEN DICYCLOPENTADIENE AND 
NITRIC ACID 


Compound C,% H, % N, % 0, % 
1 53.79 4.69 8.19 33.33 

VI 51.80 4.73 9.05 

VI 51.62 4.85 


Because of their thermal instability it was not pos- 
sible to determine the molecular weights of Solids I and 
VI by the Rast method. 

The ability of dicyclopentadiene to polymerize when 
treated with strong acids, the insolubility, and the 
amorphous character of Solids I and VI indicate that 
the latter substances are probably highly nitrated 
polymers. Dicyclopentadiene treated with concen- 
trated sulfuric acid under identical experimental condi- 
tions as with nitrie acid is readily and exothermically 
polymerized to a hard reddish-black substance. If 
several drops of sulfuric acid are added to dicyclopenta- 
diene prior to the addition of nitric acid, the ensuing 
combustion is markedly more vigorous than in the 
absence of sulfurie acid. The effect of the presence of a 
strong polymerization catalyst on the rate of the hyper- 
golic reaction can be inferred by comparing the values 
of ignition delay for various mixtures of sulfuric and 
nitric acids. Dicyclopentadiene when reacted with 
anhydrous nitric acid exhibits an ignition delay of 0.032 
sec. The addition of 5 per cent by weight of concen- 
trated (96.7 per cent) sulfuric acid to anhydrous nitric 
acid reduces the ignition delay to 0.026 sec; the addition 
of 15 per cent by weight of sulfuric acid further reduces 
the ignition delay to 0.011 sec. However, the addition 
of 5 per cent by weight of another polymerization 
catalyst, cumene hydroperoxide, served merely to dilute 
the dicyclopentadiene since, when reacted with an- 
hydrous nitric acid, the ignition delay was increased to 
0.048 sec. 


Experimental 


Preparation of Solids I and VI: Solid I was prepared 
by the dropwise addition of 0.1 mol of 99 per cent nitric 
acid to 0.1 mol of dicyclopentadiene dissolved in 60 ml 
of carbon tetrachloride and contained in a 250-ml 
beaker equipped with a mechanical stirrer. The addi- 
tion of nitric acid was regulated to keep the tempera- 
ture of the reaction between 15-20 C; if the acid is 
added too rapidly the reaction temperature rises rapidly 
and an explosive decomposition reaction may occur. 
The reddish-brown solid precipitated during the addi- 
tion of nitric acid was removed by filtration, and then 
washed with fresh portions of carbon tetrachloride and 
anhydrous ethy] ether to remove any traces of occluded 
nitric acid. The dry weight of Solid I was 3.8 g; de 
composition point 155-160 C. 
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The procedure employed for preparing Solid I was 
utilized for the preparation of Solid VI. The yield of 
solid VI was 9.6 g when 0.2 mol of nitric acid was re- 
acted with 0.1 mol of dicyclopentadiene. The decom- 
position point of Solid VI is 172 C. 

Solids I and VI do not possess melting points as 
ordinarily determined because when heated slowly they 
gradually decompose. However, there is a definite 
temperature at which Solids I and VI, when dropped 
onto a hot metallic surface, immediately decompose 
giving off a dense white vapor and oftentimes bursting 
into fame. Below the decomposition temperature, 
Solids I and VI slowly darken becoming black, but 
there is no other outward indication of their decomposi- 
tion. 

Dicyclopentadiene: The dicyclopentadiene used in all 
experiments with nitric acid was obtained as a com- 
mercial research sample. The material as received was 
a straw-colored liquid. Prior to use, the dicyclopenta- 
diene was distilled under reduced pressure to yield a 
water-white liquid boiling at 64 C (19 mm); d4”* 0.9395; 
Ip” 1.496; molecular refraction caled 40.9, found 41.1. 
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Jet Propulsion News 


By C. F. WARNER, Purdue University, Associate Editor 


Missiles 


HE Navy Department has placed an order with 

Consolidated Vultee to start producing its Terrier 
guided missile for shipboard launching against aircraft. 
It also ordered Douglas Aircraft to begin production of 
its air-to-air guided missile Sparrow. 

The Convair Terrier is a two-stage rocket developed 
by the Bureau of Ordnance. The solid-propellant mis- 
sile weighs more than a ton and a half. including the 
booster and warhead and has a range of approximately 
ten miles. The Douglas Sparrow, weighing 280 Ib 
including its warhead, was developed by the Bureau of 
Aeronautics and is designed to be used by carrier-based 
aircraft in action against enemy aircraft at ranges of 
approximately five miles. 


Aircraft 


ON July 3rd, the Navy announced that the Douglas 
Skyrocket, a supersonic research plane, has flown faster 
and higher than any other known craft. The plane was 
released from its mother ship, a B-29, over the Cali- 
fornia desert. 

Bell Aircraft Corporation has just delivered its X-5 
jet-propelled airplane to the Edwards Air Force Base, 
Muroc, Calif., where it will be flight tested. The X-5 
is the first airplane to be equipped with movable wings 
permitting the pilot to change the degree of sweepback 
during flight. For take-off, early climb, and landing, the 
wings are moved forward to obtain best performance, 
whereas the wings are moved rearward when great speed 
at high altitudes is desired. 

Northrop Aircraft, Inc., is now in quantity production 
of its all-weather interceptor, the Scorpion F-89. The 
aircraft, carrying a crew of two, is powered by two 
Allison J-35 engines equipped with afterburners. The 
F-89 has a wing span of 50 ft, is approximately 50 ft 
long, and is designed for a gross weight of over 30,000 
lb. The Air Force has announced that the Scorpion is 


NORTHROP SCORPION F-89 


BELL X-) PIONEER IN STUDY OF MOVABLE WINGS 


capable of operating at speeds in the 600-mile-per-hour 
range and of operating at altitudes in excess of 40,000 
ft. The air-borne electronic equipment enables the 
Scorpion to fly in darkness, fog, and storms, and to 
locate hostile airplanes under adverse weather condi- 
tions. 


Turbojet Engines 


THE General Electric Company has announced the 
completion of its new turbojet engine, the J-47-GE-21, 
which although the same size as its J-47 (rated in excess 
of 5200-lb thrust), is far more powerful. The new ‘‘21” 
is an all-weather engine with anti-icing features and 
high-altitude starting characteristics. It may be 
equipped for either water/alcohol injection or after- 
burning for greatly added thrust power for short 
periods. A tail-pipe variable jet nozzle may be added 
to obtain maximum efficiency. The engine has a so- 
called “cannular” type of combustor consisting of a 
single combustion space containing individual cans. 
The General Electric Company considered development 
of the compressor so significant that it bestowed its 
highest employee honor, the Coffin Award, upon two of 
the engineers responsible for the compressor’s design, 
Russell Hall and William Cornell of the Aircraft Gas 
Turbine Division. 

+ + + 

PRATT and Whitney’s newest and largest jet engine, 
the J-57, was recently test flown in a B-50. The B-50’s 
four Wasp Major engines were idled and the J-57 turbo- 
jet engine alone flew the B-50 at 370 mph. Frederick 
Brant Rentschler, president of United Aircraft Corpora- 
tion, said that Pratt and Whitney’s new J-57 has more 
power than any other turbojet engine ever flown and 
has gone a long way in reducing the enormous fuel con- 
sumption of jet engines. 


New Facilities 


A NEW plant for producing jet engine components 
will be built on the outskirts of Columbus, Ohio, by the 
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Westinghouse Electric Corporation. The plant will be 
the largest single plant yet built by the company. 


James Forrestal Research Center 


ANNOUNCEMENT was made Monday, May 28, 
1951, of the formation of a committee to raise $2,000,000 
in endowment funds for the James Forrestal Research 
Center at Princeton University, Princeton, N. J., a 
memorial to the late Secretary of Defense. 

Initial members of the committee, to be composed of 
friends and close associates of Mr. Forrestal, were an- 
nounced at a dinner at the Union Club given by Judge 
Robert P. Patterson, the committee chairman. 

In addition to the $1,500,000 provided by Princeton 
for the purchase of the property, $2,000,000 is needed 
for endowment to permit tenure of appointments to 
key scientists and for the improvement of facilities. 

The Center, devoted to co-operative enterprise in the 
fields of engineering, mathematics, chemistry, and phys- 
ies, will bring together eminent scientists and skillful 
research staffs to solve problems of prime importance 
to our national defense, transportation, and industry. 
Defense agencies have testified to their vital interest in 
the project by pledging all possible assistance and initial 
support to the program, and projected research con- 
tracts assure capacity operation of the Center’s facilities 

The James Forrestal Research Center is housed in 
two three-story buildings and 16 auxiliary buildings 
located south of U. S. Route 1, about two and a half 
miles from the center of Princeton campus. The proj- 
ects in progress, some of which are now being moved to 
the new location, are co-ordinated with the research and 
defense program of the nation. Among the projects 
planned for the new center are several important re- 
searches in jet and rocket propulsion, supersonic aero- 
dynamics, combustion, and new fuels. 


Personalities 
THE appointment of Thomas G. Lanphier, Jr., as 
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assistant to the president of Consolidated Vultee Air- 
craft Corporation, has been announced by President 
La Mott T. Cohu. Mr. Lanphier will participate in 
policy planning in the guided-missile field of the com- 
pany. Convair will build a new plant in Pomona, Calif., 
for the manufacture of guided missiles. 


The appointment of John H. Wallace, Jr., as assistant 
director of the James Forrestal Research Center at 
Princeton University, was announced by President 
Harold W. Dodds. He will assist with the development 
and administration of the laboratory facilities and with 
the organization of the joint services. The director of 
the Center is Prof. Daniel C. Sayre, formerly chairman 
of the Department of Aeronautical Engineering. 

The board of directors of Reaction Motors, Inc., 
announced that Mr. Lovell Lawrence, Jr., formerly 
president of the corporation, has been appointed chair- 
man of the board of directors. Mr. Raymond W. 
Young, former vice-president of engineering for the 
Wright Aeronautical Corporation, and later for Reac- 
tion Motors, Inc., was appointed president and general 
manager. 


American Rocket Society News 


By A. F. BOCHENEK, Associate Editor 


First ARS Meeting in Canada 
Attracts Good Attendance 


spected issues of the JouRNAL and asked 
questions about the Society’s organization, 
history, and leaders. 

Interest in the Society was also evident 


HE technical sessions and the exhibit 

of ARS literature sponsored by the 
American Rocket Society at the semi- 
‘anual meeting of The American Society 
of Mechanical Engineers held at Royal 
York Hotel, Toronto, Ontario, Canada, 
June 11-14, 1951, won many Canadian 
fiends for the Society, according to C. W. 
Chillson, Vice-President ARS, and chair- 
man of the program committee. 


SepreMBER, 1951 


Both sessions were well attended and 
in addition to prepared discussion, extem- 
poraneous comment on various aspects of 
the papers added much to the interest of 
the meetings. 


Exhibit Attracts Attention 


The exhibit of ARS literature gave an 
opportunity for many Canadians to learn 
of the work of the Society, as they in- 


from newspaper coverage of ARS sessions, 
and from radio news comment carried on 
two broadcasts. 

At the morning session, attended by 35, 
Lowell N. Randall, rocket test supervisor, 
Curtiss-Wright Corporation, Caldwell, 
N. J., discussed the operation and details 
of construction of the cavitating Venturi 
for rocket applications. Mr. Randall said 
that the cavitating Venturi has been ap- 
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plied extensively as a simple and accurate 
flow-control device, and also has been used 
for temperature measurement purposes. 

The second paper was ‘Equipment for 
Handling High-Strength Hydrogen Per- 
oxide,” by Noah S. Davis, Jr., manager, 
special products department, and John H. 
Keefe, Jr., project supervisor, Buffalo 
Electro-Chemical Company, Inc., Buffalo, 
N.Y. 

Hydrogen Peroxide Equipment 

The authors described the search made 
for the most satisfactory equipment for 
handling hydrogen peroxide in storage and 
propellant systems. Compatibility tests 
were run, they reported, on plastics, and 
various metals. They also disclosed use- 
ful information about equipment design 
and preparation of equipment for handling 
this rocket fuel. 

The final paper of the morning session 
was “Optical Methods of Rocket-Motor 
Evaluation,” by Kurt R. Stehling, rocket- 
research engineer, Bell Aircraft Corpora- 
tion, Buffalo, N. Y. 

Mr. Stehling described various methods 
tried by his company for measuring and 
analyzing thrust-chamber phenomena, 
particularly the infrared-pyrometry and 
line-reversal methods of temperature 
measurement. He also commented on 
spectrometric and  spectrophotometric 
methods as compared with conventional 
analysis methods. 


Uncooled Rocket Motors 


At the afternoon session attended by 
75, three papers were presented under the 
general subject of uncooled rocket motors. 
The Society was fortunate in having as 


chairman for this session, Lt. Col. R. A. 
Jones, AMC Power Plant Laboratory, 
Wright-Patterson Air Force Base, Dayton, 
Ohio. 


William R. Sheridan, rocket section 


engineer, Bell Aircraft Corporation, 
Buffalo, N. Y., opened the session with a 
paper on ‘‘Materials for Use in Uncooled 
Liquid-Propellant Rocket Motors.” He 
described the desirable properties for 
materials used as the refractory liners, the 
load-carrying metal outer shell, and the 
intermediate layer of insulating cement in 
construction of uncooled liquid-propellant 
rocket combustion chambers. 


Ceramic Linings 


Ceramic linings were discussed in fur- 
ther detail in the second paper, ‘“‘The 
Success of Ceramic-Lined Rocket Motors,” 
by H. Zane Schofield, supervisor, and 
W. H. Duckworth, assistant supervisor, 
ceramic division, Battelle Memorial Insti- 
tute, Columbus, Ohio. The authors re- 
ported tests which showed that sound 
practices in engineering could remove the 
observed shortcomings of such linings in 
most cases. It was their opinion that 
ceramics had a definite future in rocketry 
if their properties were properly exploited. 

In the final paper of the afternoon, 
“Silicon Carbide Linings for Uncooled 
Rocket Motors,” by Kenneth C. Nichol- 
son, supervising engineer, research and 
development department, The Carborun- 
dum Company, Niagara Falls, N. Y., 
silicon carbide as a lining material was 
discussed. Mr. Nicholson outlined the 
properties of this material and commented 
on available sizes and shapes and methods 
of fabrication and installation. 


Second International Congress on 
Astronautics Meets in London, Sept. 5—8 


THE American Rocket Society will be 
represented at the Second Inter- 
national Congress on Astronautics by 
Andrew G. Haley, general counsel ARS, 
and Lieut. Comdr. F. C. Durant III, 
members of the ARS Board of Directors. 
The congress, which will take place in 
Caxton Hall, London, England, Sept. 3-8, 
1951, will complete plans for an Inter- 
national Astronautical Federation and 
present a program of papers on the general 
theme of the Earth-Satellite Vehicle by 
international astronautical authorities. 


ARS Contribution to Program 


The American Rocket Society is spon- 
soring one technical paper. It is: “The 
Importance of Satellite Vehicles in Inter- 
planetary Flight,” by Dr. W. von Braun. 

The first part of the Congress will be 
devoted to procedural matters pertinent 
to organizing the International Astro- 
nautical Federation. These sessions will 
be open to official delegates only. Draft 


proposals for the new body have already 
been considered by the ARS Board of 
Directors and Mr. Haley, as official 
spokesman for the Society, will convey 
the ARS point of view. 

It is planned to inaugurate the new 
Federation on Tuesday, September 4. 
The next day will be devoted to detailed 
discussion for plans for the 1952 Congress 
and a program for the Federation, which 
will probably include such matters as 
standardization projects for symbols and 
nomenclature. 


16 Papers Scheduled 


The technical sessions will begin on 
Thursday morning. Sixteen papers are 
scheduled for presentation in four sessions. 
While the public will be invited, the pro- 
gram is expected to have little interest for 
laymen. 

The program will conclude with a 
popular technical symposium on Saturday 
afternoon at which spokesmen from the 


various national delegations will give 
short talks explaining their ideas on inter. 
planetary flight. A public statement wil] 
also be made giving the results of the 
preceding Congress sessions. 

Some of the other papers on the tech- 
nical program are: Applications of Satel- 
lite Vehicles, by E. Birch Andersen (Den- 
mark); The Ascent of Satellite Vehicles, 
by D. F. Lawden (Britain); Views on the 
Ascent of the Earth Satellite Vehicle, by 
Ing. Janson (Germany); Start, Return, 
and Landing of an Optimum Satellite Step 
Rocket, by Dipl. Ing. H. Kuhme (Ger. 
many). 

As the final event of the Congress, 
delegates will be entertained at dinner at 
the St. Ermin’s Hotel, tendered by the 
British Interplanetary Society. 


Student Papers Heard at 
New York Section 
Meeting 


NDERGRADUATE research _ proj- 
ects carried out by students at Pratt 
Institute, Brooklyn, N. Y., and at the Col- 
lege of the City of New York, New York, 
N. Y., were described at a meeting of the 
New York Section of the American Rocket 
Society held at the Engineering Societies 
Building, New York, N. Y., May 25, 
1951. The meeting was arranged by 
Major James R. Randolph, a member of 
the Board of Directors of the New York 
Chapter. Fifty members and guests were 
present. Major Randolph presided. 

Design and construction of an improved 
nozzle and test section and the use of 
superheaters to obtain a high Mach 
number flow in a supersonic wind tunnel 
were described by Roy O. Bohner and 
Vincent P. Buschimi, Pratt Institute 
students. 

The second talk was by V. A. Chobotov 
and Howard Dobias of Pratt Institute. 
They described the design of a pulse jet 
to be used for static testing of fuels, 
nozzles, and other factors involved in 
pulse jets. 

Ralph A. Benson of Pratt Institute de- 
scribed his work on a static test rocket for 
gaseous fuels. He pointed out that his 
rocket could be used for testing propellants 
and mixing devices. It was fired upward, 
was water cooled, and was designed to be 
used with an ordinary household scale to 
measure thrust. 

The last project described was by 
Howard Macklis, Sheldon Newberger, and 
Harris L. Shapiro of the College of the City 
of New York. They demonstrated how 
the theory of gas turbines is applied to 4 
specific design problem. Their design gave 
a thrust of 2600 Ib. They showed also 
how dimensions and angles of blades and 
other details not derived from thermo 
dynamics were calculated. 
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Indiana Section Completes Memorable 
First Year 


HE Indiana Section of the American 

Rocket Society climaxed a memorable 
first vear of operation by an “Acceleration 
Banquet”’ held in the North Ballroom of 
the Purdue Memorial Building, West 
Lafayette, Ind., May 24, 1951. A dis- 
tinguished audience of nearly 100 persons 
which included officials of the University, 
Lt. Gov. John A. Watkins of Indiana, and 
H.R. J. Grosch, president of the American 
Rocket Society, heard Under-Secretary of 
the Navy Dan A. Kimball talk on the history 
and present status-of rocket development. 
(Mr. Kimball has since been named 
Secretary of the Navy.) J. P. Layton, 
president of the Indiana Section, presided. 

Mr. Layton introduced charter members 
of the Section, newly elected officers, and 
members of the Board of Directors of the 
Section and reviewed the program of the 
past season. 

Purdue Officials Speak 

Following a brief introduction by Presi- 
dent Layton, R. B. Stewart, vice-presi- 
dentof Purdue University, who represented 
Purdue’s president, F. L. Hovde, expressed 
the good wishes of the University for the 
continued success of the Section and con- 
gratulated the Section for its recognition 
by the faculty as a university activity. 
He was followed by the Honorable John 
A. Watkins, lieutenant-governor of In- 
diana, who represented Governor Schricker. 
Lt. Gov. Watkins spoke highly of the 
scientific work of the ARS. He was fol- 
lowed by A. A. Potter, dean of the Univer- 
sity’s School of Engineering, who reviewed 
the development of engineering education 
in the United States, particularly the con- 
tributions of the land-grant colleges. 
Dean Potter cited many “‘firsts’”’ in engi- 
neering education at Purdue, and then 
stressed the importance of the humanities 
as well as technology for the development 
of good character and responsible citizen- 
ship. 

ARS President Present 

President Grosch spoke briefly to re- 
view the development of the American 
Rocket Society from its early days of back- 
yard experimentation to its present status 
is one of the leading technical organiza- 
tions working in the field of rocket and 
iet-propulsion technology. Rocketry, he 
commented, was a young man’s field and 
one in which talented hard-working young 
men can progress far in a few years. 
President Grosch announced expansion of 
the ARS Journat both in format and in 
lumber of issues. He spoke appreciatively 
of the Society’s joint activities with The 
American Society of Mechanical Engineers 
ind the Institute of the Aeronautical 
kiences. 
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SECRETARY OF THE NAVY DAN A. KIMBALL, 
MAIN SPEAKER AT THE ARS INDIANA SEC- 
TION DINNER HELD May 24, 1951 


The Society’s Board of Directors, he 
said, was working to increase aid on the 
national level to the ARS Sections and to 
organize Sections in centers of rocket re 
search and industry not now being served 
by Sections. The ARS Fall meeting at 
Minneapolis, Minn., he added, would give 
members in the Midwest an opportunity 
to attend a major meeting of the Society. 
He invited all to attend and to participate 
in the discussion of the program of inter- 
esting papers. 


Promising Progress 


Following a welcome from Capt. J. C. 
Woelfel of the Purdue Naval Reserve 
Officers Training Corps, President Layton 
introduced the Honorable Dan A. Kimball, 
Undersecretary of the Navy. In spite of 
promising progress in the field of rockets 
and jet propulsion, Mr. Kimball said, a 
great deal remains to be done. Ever- 
increasing technological challenges face 
members of the ARS. As development 
progresses, rockets will become more and 
more useful instruments of national de- 
fense and civil progress, he predicted. 


New Officers Elected 


At a meeting of the Indiana Section on 
May 3, 1951, the following officers were 
elected: President, David MeNay; 
Vice-President, Clair M. Beighley; Secre- 
tary, Vincent Capasso; Treasurer, Ralph 


Kester; Board of Directors, M. J. Zucrow, 
W. Gilliland, Arthur Wiggins, Gladys 
Geiger, and J. P. Layton. 


ARS Southern California 
Section Holds Closed 
Meeting 


following officers have been elected 

for the 1951-1952 season by the 
Southern California Section of the Ameri- 
can Rocket Society: President, R. B. 
Canright, Jet Propulsion Laboratory, 
California Institute of Technology; Vice- 
President, B. L. Dorman, Aerojet Engi- 
neering Corporation; Secretary-Treasurer, 
R. J. Lodge, North American Aviation, 
Inc. Directors elected were as follows: 
C. C. Ross, Aerojet Engineering Corp.: 
L. G. Dunn, Jet Propulsion Laboratory, 
California Institute of Technology; J. F. 
Manildi, University of California, Los 
Angeles, Calif.; G. D. Brewer, Hughes 
Aircraft Corp.; W. T. Cox, USN, Bureau 
of Aeronautics; W. E. Zisch, Aerojet 
Engineering Corp.; E. G. Crofut, Aerolab 
Development Company. 

The Southern California Section held a 
highly successful closed dinner meeting at 
the Pasadena Athletic Club on April 24, 
1951. The title of the meeting was 
“Symposium on Propellant Selection.’ 
President Canright led the discussion by a 
panel of rocket authorities which included: 
Tom Dixon, North American Aviation, 
Inc.; Don L. Armstrong, Aerojet Engi- 
neering Corp.; Leland G. Cole, Jet Propul- 
sion Laboratory, California Institute of 
Technology; David A. Young, Aerojet 
Engineering Corp.; and T. E. Myers. 
North American Aviation, Inc. 

The meeting was of such interest that 
few of the nearly 200 in the audience left 
before 11:30 p.m. The meeting was also 
distinctive because it was the largest 
classified meeting ever held by the Section. 
In addition to local members and guests, 
all local military services were well repre- 
sented. The meeting closed with the 
showing of a film from North American 
Aviation, Inc., on testing propellants in a 
small engine. 


Interesting Papers 
Planned for 1951 
ARS Fall Meeting 


HE American Rocket Society will hold 

its 1951 Fall Meeting in conjunction 
with The American Society of Mechanical 
Engineers at the Radisson Hotel, Minne- 
apolis, Minn., Sept. 25-28, 1951. 

The ARS part of the program will be a 
session on Thursday afternoon, September 
25, during which two papers will be pre- 
sented. The first will be by Clark E. 
Thorp, chairman of chemistry and chemi- 
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cal engineering research, Armour Re- 
search Foundation, Illinois Institute of 
Technology, Chicago, Ill. Dr. Thorp’s 
subject will be “The Properties and Pro- 
duction of Concentrated Gaseous and 
Liquid Ozone.” 

The second paper by Stanley V. Gunn, 
Purdue University Rocket Laboratory, 
Lafayette, Ind., will be on “The Effects of 
Several Variables Upon the Ignition of 
Hypergolic Fuels Oxidized by Nitric 
Acid.”” Mr. Gunn has divided his paper 
into three parts covering (1) description 
of electronic timer and reaction apparatus 
employed for open-cup tests; (2) igni- 
tion-lag data obtained from open-cup 
tests; and (3) a description of the reac- 
tion apparatus to be employed for igni- 
tion-lag determinations within a small 
rocket motor. He will also report data 
which show the effects of propellant tem- 
perature and water content of the oxidizer 
on ignition lag of several hypergolic fuels, 
when oxidized with nitric acid. 

Members of the Society, especially 
those residing in the Midwest, are urged 
to participate in the Fall Meeting. Engi- 
neers who are not ARS members but who 
are working in the rocket industry are also 
invited to the ARS session. 


New York Section Plans 
Fall Meetings 


IHE ARS New York Section will open 

its 1951-1952 series of meetings at the 
Engineering Societies Building, New York, 
N. Y., Sept. 21, 1951, when Edward F. 
Chandler, engineer and inventor, of Brook- 
lyn, N. Y., will talk on “Adapting the 
Rocket to Modern Warfare.”’ 

Mr. Chandler will review some phases 

of the work he is doing on a method of 
firing rocket-type ammunition from 
breech-loading firearms. He is of the 
opinion that the tactical uses of the rocket 
can be increased considerably by designing 
a rocket as a self-contained round or 
“fixed ammunition.” 
* This would make available light, highly 
mobile arms capable of being rapidly 
served and fired either manually or auto- 
matically. 

The regular meeting will be held on 
Friday, Oct. 19, 1951, at 8 p.m. The 
speaker will be Dr. R. W. Porter, Divi- 
sion Engineer, Special Project Division, 
General Electric Company, Schenec- 
tady, N. Y. The subject of Dr. Porter’s 
talk will be ‘Failures in the V-2 Pro- 


Dr. Porter will discuss the V-2 
firing with particular regard to material 
difficulties which prevented accomplish- 
ment of specific objective desired. 


gram.” 


ARS Junior Award 


Competition 


(He distinction of receiving an Ameri- 

can Rocket Society Junior Award car- 
ries with it a recognition of meritorious 
achievement early in one’s professional 
career. Papers are now being accepted 
at the New York office for consideration 
for the 1951 award. The award and 
medal will be presented at the Annual 
Convention to be held in November. To 
be considered, papers must be received 
not later than Oct. 20, 1951. 

The winning paper will be judged mainly 
on the basis of content which should re- 
flect original thought and effort. The 
age of authors of papers must not be 
greater than 25 years. Papers should be 
on standard size paper, typewritten, and 
clearly marked, “Submitted for Junior 
Award Competition.’’ Send papers to: 
Secretary, American Rocket Society, 29 
West 39th Street, New York 18, N. Y. 


Structures Engineers 


Dynamicists 


WRITE 


P.O. Box 516 


ENGINEERS 


Supersonic Aerodynamicists 


Heat Transfer Specialists 
Flutter and Vibration Specialists 
Servomechanism Designers 


for Missile Design and 


Development Work 


TECHNICAL PLACEMENT OFFICE 
McDONNELL AIRCRAFT CORPORATION 


ST. LOUIS 3, MISSOURI 


SELF-ACTUATED 


CONTROLS AND REGULATORS. 


High pressure liquids and gases 


Pressure Reducing 
Back Pressure 


Special light weight units 
for air-borne service. 


Special materials for corrosive 
or low-temperature applications. 


GROVE REGULATOR COMPANY 


Oakland 


Los Angeles 


GROVE 


for 


Relief and By-Pass 
Manual Valves 


Houston 


New York 


ARS Journal 


Pp. 


el 
conce 
and ec 
teen 
back 
comple 
form. 

‘spects 
thirty - 
for con 


22 dise 
hop co 
the cor 
ind th 
their a 
hop 
thapter 
¥rvome 
thapter. 
There 
he end 

7 
‘aching 
thors 
lect rie 
iantita 


DEPT EM 


E 
no 
ra 
to 
rea 
3 bo 
of 
trit 
whi 
well 
ma) 
lead 
mar 
Ir 
issue 
some 
year 
issue 
Rol 
|| 
GROVE | 
SIMPLICITY 
ACCURACY 
SAFETY | 
| 


sh- 


Technical Literature Digest 


By H. S. SEIFERT, California Institute of Technology, Associate Editor 


CONTRIBUTORS: 


F. C. Gunther 
D. Altman 
S. S. Penner 


R. B. Canright 
D. I. Baker 
H. S. Siefert 


Epiror’s Note: 

The following collection of references is 
not intended to be comprehensive, but is 
rather a selection of the most significant 
and stimulating papers which have come 
to the attention of the contributors. The 
reader will understand that a considerable 
body of literature is unavailable because 
of security restrictions. We invite con- 
tributions to this department of references 
which have not come to our attention, as 
well as comment on how the department 
may better serve its function of providing 
leads to the jet-propulsion applications of 
many diverse fields of knowledge. 

Inasmuch as this represents the first 
issue of the Literature Digest department, 
wme references dated even more than one 
year ago have been listed. In future 
issues, the contributors plan to devote the 
space primarily to current listings. 


Book Review 


SERVOMECHANISMS AND REGULATING Sys- 
Tt—M Desicn, by Harold Chestnut and 
Robert W. Mayer, John Wiley & Sons, 
Inc., New York, N. Y., vol. I, xiii + 505 
pp. $7.75. 


HIS book is a carefully executed 

engineering text on a subject of vital 
concern to those interested in the guidance 
ind control of missiles. Three of its four- 
ten chapters review useful mathematical 
wwekground, including applications of 
complex numbers and the Laplace trans- 
orm. There is a brief review of certain 
‘spects of network theory, followed by a 
thirty-page discussion of stability criteria 
ir control systems. Chapters 7 through 
12 discuss many examples of single closed- 
hop control systems (including on p. 190 
ie control of a missile in vertical flight) 
ind the use of graphical techniques for 


heir analysis and synthesis. Multiple- 
bop control systems are discussed in 
thapter 13, and the correlation between 
lansient and steady-state performance of 
*tvomechanisms is examined in the last 
thapter. 
There is a fifty-page set of problems at 
le end of the book, a very valuable fea- 
This book is based on extensive 
“aching and engineering experience of the 
uthors and their colleagues at the General 
Hketric Co., and is a thorough, sound, and 
antitative approach to an important 
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subject. Volume II of the set will extend 
the subject to more detailed consideration 
of practical design problems, synthesis of 
systems, “noise,” nonlinearities, and the 
like. 


Books 


The Behavior of Engineering Metals, 
by H. W. Gillett, John Wiley & Sons, 
Inc., New York, N. Y., 1951. 

Engineering Supersonic Aerodynamics, 
by E. A. Bonney, McGraw-Hill Book Co., 
Inc., New York, N. Y., 1950 (pp. 60-69, 
246-252—Introduction to ramjet prin- 
ciples). 

A World Airlift, by Elvira K. Frodkin, 
Funk & Wagnalls Co., New York, N. Y., 
1950, viii + 216 pp. $2.85. 

Heat and Temperature Measurement, 
by R. L. Weber, Prentice-Hall, Inc., New 
York, N. Y., 1950, 422 pp. 

Cloud Physics, by D. W. Perrie, John 
Wiley & Sons, Inc., New York, N. Y., 
1950, 30 plates and 119 pp. $4.50. 

Ballistic Disturbance Theory (Ballis- 
tische Storungs-Theorie), by Raymund 
Sanger, Verlag Birkhauser, Basel, Switzer- 
land, 1949, xii + 226 pp. Swiss francs, 
$14.50. 

Numerical Analysis of Heat Flow, by 
G. M. Dusinberre, McGraw-Hill Book Co., 
Inc., New York, N. Y., 1949, xiv + 227 pp. 

Heat Transfer, Vol. I, by Max Jakob, 
John Wiley & Sons, Inc., New York, N. Y., 
1949. 

Helicopter Engineering, by Raymond 
A. Young, The Ronald Press Co., New 
York, N. Y., 1949, xiv + 255 pp. $10 
(pp. 130-134, Discussion of jet-propelled 
helicopter rotors). 

Rockets, Missiles, and Space Travel, 
by Willy Ley, Viking Press, New York, 
N. Y., 1951. 


Jet Propulsion Engines 


Burners for Supersonic Ramjets, by 
J. W. Mullen, J. B. Fenn, and R. C. 
Garmon, Industrial and Engineering Chem- 
istry, vol. 43, no. 1, Jan., 1951, p. 195. 

Method for Calculation of Ramjet 
Performances, by J. R. Henry and J. 
Buel Bennett, NACA Technical Note 
N2357, June, 1951. 

Flame Blowoff Studies of Cylindrical 
Flame Holders, by G. W. Haddock, 
Jet Propulsion Laboratory, Pasadena, 
Calif., PR 3-24, May, 1951. 


High Speed Aerodynamic Problems of 
Turbojet Installations, by H. Luskin and 
H. Klein, Douglas Aircraft Co., Santa 
Monica, Calif., Report 4SM-13830 (pre- 
sented at the ASME Annual Meeting, 
Dec., 1950, published Sept., 1950). 

Gas Dynamic Theory of the Ramjet, 
by J. Reid, Aeronautical Research Coun- 
cil, London, Report Memorandum 2370, 
1950, p. 52. 

Combustion Problems in Ramjet Design, 
by J. P. Longwell, Journal of the Aero- 
nautical Sciences, vol. 16, no. 12, Dec., 
1949, pp. 707-714. 


Rocket Propulsion 
Engines 


New Techniques for Obtaining Heat 
Transfer Parameters of the Wall and Com- 
bustion Gas in a Rocket Motor, by M. E. 
Ellion, Transactions of The American 
Society of Mechanical Engineers, vol. 73, 
no. 2, Feb., 1951, pp. 109-114. 

Influence of After-Burning on the Per- 
formance of Rockets (in French), by R. 
Engel and U. T. Boedewadt, La Recherche 
Aéronautique, no. 18, Nov.-Dec., 1950, pp. 
23-33. 

Stability of Flow in a Rocket Motor, 
by D. F. Gunder and D. R. Friant, 
Journal of Applied Mechanics, vol. 17, 
Sept., 1950, pp. 327-333. 

Early History of Rocket Research, by 
J. W. Siry, Scientific Monthly, vol. 71, 
Nov., 1950, pp. 326-332. 

Rocket Research in the Twentieth 
Century, by J. W. Siry, Scientific Monthly, 
vol. 71, Dec., 1950, pp. 408-421. 

Cooling of the Rocket Combustion 
Chamber (in German), by H. Gartmann, 
Weltraumfahrt, vol. 5, Oct., 1950, pp. 109- 
114. 

Is the Nuclear-Powered Rocket Feas- 
ible? by H. S. Seifert, Physics Today, 
vol. 2, May, 1949, pp. 6-13. 

Preliminary Experiments of Gaseous 
Transpiration Cooling of Rocket Motors, 
by R. B. Canright, Jet Propulsion Labora- 
tory, Pasadena, Calif., PR 1-75, Nov. 24, 
1948. 

Evaluating Rocket Propellant Perform- 
ance, by Joseph Himpon, /nteravia, vol. V, 
no. 10, Oct., 1950, pp. 530-532. 


Heat Transfer and 
Fluid Flow 


Heat Transfer, Pressure Drop and Burn- 
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out Studies with and without Surface 
Boiling for De-aereated and Gassed Water 
at Elevated Pressures in a Forced Flow 
System, by H. Buchberg, F. Romie, R. 
Lipkis, and M. Greenfield.! 

On the Calculation of the Rate of Heat 
Transfer Through a Laminar Boundary 
Layer, by Arthur N. Tifford and John 
Wolansky, Journal of the Aeronautical 
Sciences, vol. 18, no. 6, June, 1951. 

Heat Conduction in a Compressible 
Fluid, by J. D. Cole and T. Y. Wu.! 

Heat Transfer and Pressure Drop Data 
for High Heat Flux Densities to Water at 
High Sub-Critical Pressures, by Warren 
M. Rohsenow and John A. Clark.! 

Internal Flow Data and a Heat Transfer 
Theory for the Vortex Refrigerating Tube, 
by George W. Scheper, Jr.! 

The Transient Method for Determining 
Heat Transfer Conductance from Bodies 
in High Velocity Fluid Flow, by C. R. 
Garbett.! 

Heat Transfer in Turbulent Shear 
Flow, by W. D. Rannie, thesis, Guggen- 
heim Aeronautical Laboratory, California 
Institute of Technology, Pasadena, Calif., 
June, 1951. 

Variation of the Eddy Conductivity 
With Prandtl Modulus and Its Use in 
Prediction of Turbulent Heat Transfer 
Coefficients, by Rodman Jenkins.! 

Liquid Superheat and Boiling Heat 
Transfer, by B. R. Mead, F. E. Romie, 
and A. G. Guibert.! 

Some Experiments on the Heat Transfer 
from a Gas Flowing Through a Conver- 
gent-Divergent Nozzle, by O. A. Saunders 
and P. H. Caulder.! 

Photographic Study of Surface-Boiling 
Heat Transfer to Water with Forced 
Convection, by F. C. Gunther, T’ransac- 
tions of The American Society of Mechani- 
cal Engineers, vol. 73, no. 2, Feb., 
1951, pp. 115-123. 

Approximate Solutions for Shock Waves 
in a Steady, One-Dimensional, Viscous 
and Compressible Gas, by P. Lieber, F. 
Romano, and H. Lew, Journal of the 
Aeronautical Sciences, vol. 18, no. 1, Jan., 
1951, pp. 55-60. 


Combustion 


On the Theory of Combustion of 
Initially Unmixed Gases, by Y. B. Zeldo- 
vich, NACA Technical Memorandum 
1296, June, 1951. 

Influence of Flame Front on the Flow 
Field, by H. 8. Tsien, Journal of Applied 
Mechanics, vol. 18, June, 1951, p. 188. 

Theory of Flame Propagation, by Y. B. 
Zeldovich, NACA Technical Memorandum 
1282, June, 1951. 

Interaction of Flow Pulsations and 
Flame Propagation, by G. H. Markstein, 
Journal of the Aeronautical Sciences, vol. 
18, 1951, p. 428. 

Flames: Their Structure and Radiation, 
by A. G. Gaydon, Endeavour, vol. 10, 
Jan., 1951, pp. 17-21. 


1 1951 Heat Transfer and Fluid Me- 
chanics Institute (preprints of papers pre- 
sented at Stanford University, Palo Alto, 
Calif.), June 20-22, 1951. 
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On the Burning Velocity of Double-Base 
Powders in Nitrogen Under Pressures 
From 100 to 10,000 Kg/Cm? (in French), 
by H. Muraour and J. Faureau, Chemie et 
Industrie, vol. 65, no. 1, Jan., 1951, pp. 
53-55. 

On the One-Dimensional Theory of 
Flame Structure, by R. Friedman and E. 
Burke, Journal of the Aeronautical Sciences, 
vol. 18, no. 4, April, 1951, pp. 239-246. 


Preliminary Studies on the Ignition and 
Propagation of Flames in High-Speed 
Gas Streams, by H. L. Olsen and E. L. 
Gayhart, The Johns Hopkins University, 
Silver Spring, Md., Report No. CM-632, 
Feb. 8, 1951. 

Flame Propagation II. The Determina- 
tion of Fundamental Burning Velocities of 
Hydrocarbons by a Revised Tube Method, 
by M. Gerstein, O. Levine, and E. L. 
Wienz, Journal of the American Chemical 
Society, vol. 73, Jan., 1951, pp. 418-422. 


Flame Propagation III: Theoretical 
Consideration of the Burning Velocities of 
Hydrocarbons, by D. M. Simon, Journal 
of the American Chemical Society, vol. 73, 
Jan., 1951, pp. 422-425. 


Symposium on Kinetics of Propellants 
(Burning Rate of Double-Base Powders), 
by R. E. Gibson, B. L. Crawford, Jr., 
et al, F. Daniels, et al, J. H. Frazer, et al, 
B. Lewis, et al, O. K. Rice, et al, W. H. 
Avery, Journal of Physical and Colloid 
Chemistry, vol. 54, 1950, pp. 847-954. 


Investigation of Turbulent Flames, by 
B. Karlovitz, D. W. Denniston, Jr., and 
F. E. Wells, Journal of Chemical Physics, 
vol. 19, 1951, p. 541. 


Flame Propagation IV. Correlation of 
Maximum Fundamental Flame Velocity 
with Hydrocarbon Structure, by R. R. 
Hibbard and B. Pinkel, Journal of the 
American Chemical Society, vol. 73, April, 
1951, pp. 1622-1625. 


Experimental and Theoretical Studies 
of Flame-Front Stability, by G. H. Mark- 
stein, Journal of the Aeronautical Sciences, 
vol. 18, no. 3, March, 1951, pp. 199-209. 


An Analog Computer for Flame Gas 
Composition, by W. 8S. McEwan and S§. 
Skolnik, Review of Scientific Instruments, 
March, 1951, pp. 125-132. 


Determination of the Average Tempera- 
turein a Combustion Chamber (in French), 
by L. Viaud, La Recherche Aéronautique, 
no. 18, Nov.-Dec., 1950, pp. 55-58. 


Theory of Combustion Limits and Igni- 
tion of Combustible Mixtures (in German), 
by K. Bechert, Annalen der Physik, vol. 7, 
no. 3-4, 1950, pp. 113-128. 


Ionization and Luminescence in Flames 
(in German), by E. Sanger, I. Bredt, and 
P. Goercke, Zeitschrift fiir angewandte 
Mathematik und Physik, vol. 1, 1950, 
pp. 209-210 (Translated: NACA Technical 
Memorandum 1305, April 1951). 


Hypothesis of Powder Burning (in 
German), by E. Sanger, Zeitschrift fiir 
Naturforschung, vol. 5a, no. 8, 1950, pp. 
467-468. 


Basic Studies on Flame Stabilization, by 
G. C. Williams, Journal of the Aeronautical 
Sciences, vol. 16, no. 12, Dec., 1949, pp. 
714-722. 

The Relationship Between Performance 


and Constitution of Pure Organic Fx. 
plosive Compounds, by W. C. Lothrop and 
G. R. Handrick, Chemical Reviews, vol, 
44, no. 3, June, 1949, pp. 419-445. 

The Theory of Flame Propagation, by 
J. O. Hirschfelder and C. F. Curtiss, 
Journal of Chemical Physics, vol. 17, no, 
11, Nov., 1949, pp. 1076-1081. 


Fuels, Propellants, and § - 
Materials 


Thermal Decomposition of Hyd) azine 
Vapor in Silica Vessel, by T. J. Hanvatty, 
J. N. Pattison, J. W. Clegg, and A. W. 
Lemmon, Jr., Industrial and Engin: ering 
Chemistry, vol. 43, May, 1951, pp. 7 113- 
1116. 

The Mechanism of the Catalytic De- 
composition of Hydrogen Peroxid> by 
Silver, by R. L. Wentworth, Mass:ichu- 
setts Institute of Technology, Cambvidge, 
Mass., Report 32, May 1, 1951. 

Methy! Alcohol: The Entropy, Heat 
Capacity and Polymerization Equilibria 
in the Vapor, and Potential Barrier to 
Internal Rotation, by W. Weltner, Jr., 
and K. Pitzer, Journal of the rican 
Chemical Society, vol. 73, 1951, p. 2606. 

Oxygen and Nitrogen Isotope Effects 
in the Decomposition of Ammonium Ni- 
trate, by L. Friedman and J. Bigeleisen, DE 
Journal of Chemical Physics, vol. 18, 5 0 
1951, p. 1325. 

Calculated Thermochemical Properties 
for Hydrazine-Hydrogen Peroxide Reac- 
tion at High Pressures, by W. J. Barr, and 
E. J. Wilson, Jr., Journal of Chemical 
Physics, vol. 19, 1951, p. 470. 

New Zirconia Refractory Material Use- 
ful at Temperatures to 4600°F., by 
Anonymous, Materials and Methods, 
March, 1951, p. 81. 

Heat Resistant Ceramic Coatings, by 
D. G. Gennett, Materials and Methods, 
March, 1951, pp. 65-67. 

The Mechanism and Kinetics of the 
Water-Gas Reaction in Cooling Combus- 
tion Gases (in German), by H. Behrens, 
Zeitschrift fiir Physikalische Chemie, vol. 
195, 1950, p. 24. 

Rocket Propellants, by Heinz Gartmann, 
Interavia, vol. V, no. 10, Oct., 1950, pp. 


527-529. 
Mea 

Means 

Physical-Chemical Topics 


Clas: 
Time Lag in the Self-Nucleation of fio R. 
Supersaturated Vapor, by R. F. Prob fManes, 
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Means of the Critical Flow Orifice, by 
i. Wise, and D. Altman, Journal of 
Chemical Physics, vol. 19, 1951, p. 515. 
Classical Thermodynamics and Reac- 
ion Rates Close to Equilibrium, by M. 
Manes, L. J. E. Hofer, and S. Weller, 
mrnal of Chemical Physics, vol. 18, 1951, 
1355. 

Surface Tension of Small Droplets as a 
function of Size from Critical Supersatura- 
im Data, by G. M. Pound and V. K. 
laMer, Journal of Chemical Physics, vol. 
N,1951, p. 506. 

High Intensity Source for the Molecular 
vam. Part I. Theoretical, by A. 
Mantrowitz and J. Grey. Part II. Ex- 
imental, by G. B. Kistiakowsky and 
.P. Slichter, Review of Scientific Instru- 
8, vol. 22, 1951, p. 328. 


Mechanism and Kinetics of the Water- 
Gas Reaction During Adiabatic Expansion 
(in German), by H. Behrens, Zeitschrift 
fiir Physikalische Chemie, vol. 195, 1950, 
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Tables of Contributions to Thermo- 
dynamic Properties due to Gas Imperfec- 
tion, by W. B. White, The Rand Corpora- 
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Thickness of Shock Fronts in Argon and 
Nitrogen and Rotational Heat Capacity 
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A Study of Slow Self-Propagating Re- 
actions in Solids, by R. A. W. Hill, L. E. 
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Research, vol. 3, 1950, p. 569. 
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272-280. 
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An Analysis of a Compound Pendulum 
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thesis, Guggenheim Laboratory, Cali- 
fornia Institute of Technology, Pasadena, 
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Instruments, by G. P. L. Williams, 
Society of Instrument Technology, Transac- 
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Use of Resistance Thermometers for 
Measuring Rapidly Changing Tempera- 
tures, by C. B. Daish, et al, Philosophical 
Magazine, vol. 41, July, 1950, pp. 729-730. 


Vortex Thermometer for Measuring 
True Air Temperatures and True Air 
Speeds in Flight, by B. Vonnegut, Review 
of Scientific Instruments, vol. 21, Feb., 
1950, pp. 136-141. 

Response of Thermocouple to Rapid 
Gas-Temperature Changes, by M. W. 
Carbon, et al, Transactions of The American 
Society of Mechanical Engineers, vol. 72, 
July, 1950, pp. 655-657. 

Sonic-Flow-Orifice Temperature Probe 
for High-Gas-Temperature Measurements, 
by P. L. Blackshear, Jr.. NACA Tech- 
nical Note 2167, Sept., 1950. 


A Simple Pressure Balance, by C. A. 
Swenson, Review of Scientific Instruments, 
vol. 21, Jan., 1950, pp. 22-24. 


Dynamic Accuracy in Temperature 
Measurement, by J. G. Ziegler and N. B. 
Nichols, nstruments, vol. 23, Jan., 1950, 
pp. 66-69. 

Heat-Conduction Errors in Temperature 
Measurement, by L. E.Smith, Transactions 
of The American Society of Mechanical Engi- 
neers, vol. 72, Jan., 1950, pp. 71-76. 

Attenuation of Oscillatory Pressures in 
Instrument Lines, by A. S._ Iberall, 
Transactions of The American Society of 
Mechanical Engineers, vol. 72, July, 1950, 
pp. 689-695. 

A Versatile Pneumatic Instrument 
Based on Critical Flow, by W. A. Wild- 
hack, Review of Scientific Instruments, vol. 
21, Jan., 1950, pp. 25-30. 

Measurement of Short Time Intervals, 
by K. J. Brimley, Transactions of the 
Society of Instrument Technology, vol. 2, 
March, 1950, pp. 2-17. 

Resistance-Wire Strain Gauges, by J. 
G. Strong, Journal of the Royal Aeronauti- 
cal Society, vol. 54, Jan., 1950, pp. 19-26. 


Flight, Ballistics, and 
Vehicle Development 
The Use of Auxiliary Rockets in High- 


Speed Aircraft, by Z. Plaskowski, Aircraft 
Engineering, March, 1951, pp. 72-75. 
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A Theoretical Analysis of the Effects of 
Fuel Motion on Airplane Dynamics, by 
Albert A. Schy, NACA Technical Note 
2280, Jan., 1951, 53 pp. 

Rockets and Assisted Take-Off, by A. V. 
Cleaver, Journal of the Royal Aeronautical 
Society, Feb., 1951, pp. 87-109. 

The Perturbation Calculus in Missile 
Ballistics, by R. Drenick, Journal of the 
Franklin Institute, April, 1951, pp. 423-436. 


Phenomena Associated with the Flight 
of Ultra-Speed Pellets. Part I: Ballistics. 
Part II: Spectral Character of Luminos- 
ity. Part III: General Character of 
Luminosity, by J. S. Rinehart, W. A. 
Allen, and W. C. White, Bulletin of the 
American Physical Society, vol. 26 (9), 
July 25-28, 1951. 

Jet Propulsion and Guided Missiles, 
by Anonymous, Anti-Aircraft Journal, 
March-April, 1951, pp. 46-53. 

The Nyquist Criterion of Stability, by 
J.C. West, Electronic Engineering, vol. 22, 
May, 1950, pp. 169-172. 

Aircraft Radio Altimeters, by B. A. 
Sharpe, Journal of the Institute of Naviga- 
tion, vol. 3, Jan., 1950, pp. 79-89. 

An Analysis of Base Pressure at Super- 
sonic Velocities and Comparison with 
Experiment, by Dean R. Chapman, 
NACA Technical Note 2137, July, 1950, 
69 pp. 

The Shape of Ships to Come, by W. Ley, 
Interavia, vol. V, no. 10, Oct., 1950, pp. 
496-499. 

Project Reach, by C. DeVore, Naviga- 


tion, vol. 2, Dec., 1950, pp. 275-281. 


Rocket Navigation, by 8S. Herrick, 
Navigation, vol. 2, Dec., 1950, pp. 259- 
272. 

Earth Satellite Vehicles, by R. Engel, 
Interavia, vol. V, no. 10, Oct., 1950, pp. 
500-502. 

Limitations of Rocket Propulsion and 
Their Significance for the Idea of Space 
Travel (in German), by Rudolf H. Reichel, 
Zeitschrift des Vereines deutscher Ingenieure, 
vol. 92 (32), Nov. 1950, pp. 873-882. 

A Speed Limit for Rockets, by W. 
Squire, Indian Journal of Physics, vol. 24, 
Sept., 1950, p. 410. 

Base Pressure at Supersonic Velocities, 
by A. Gabeaud, Journal of the Aeronautical 
Sciences, vol. 17, Aug., 1950, pp. 525-526. 

The Sperry Zero Reader, by 8S. Kellogg 
and C. F. Fragola, Aeronautical Engineer- 
ing Review, vol. 8, Nov., 1949, pp. 22-31. 

Plausible Speed Limit for Metallic 
Projectiles, by K. R. Dixit, Indian Journal 
of Physics, vol. 23, Nov., 1949, pp. 469- 
476. 

A Generalization of Nyquist’s Stability 
Criterion, by A. Vazsonyi, Journal of 
Applied Physics, vol. 20, Sept., 1949, pp. 
863-867. 
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